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1  4  A 


Conceptually,  a  barrier  weapon  system  consisting  of  a 
warhead  that  could  deliver  a  large  number  of  projectiles  at 
very  high  velocities  would  be  an  extremely  effective  method  of 
defeating  armored  vehicles  such  as  tanks  and  personnel  carriers. 
Assuming  that  the  projectiles  could  be  adequately  dispersed,  the 
principal  advantage  of  such  a  warhead  would  be  the  larger 
target  area  that  could  be  covered  with  one  warhead.  Since  the 
energy  required  to  launch  projectiles  to  a  given  velocity  is  a 
direct  function  of  the  projectile  mass,  twenty  times  as  many 
5-gram  projectiles  could  be  launched  to  a  specific  velocity  as 
a  single  100-gram  projectile.  Thus,  for  a  given  launcher  weight 
and  volume,  a  much  larger  area  could  be  controlled  by  the  use 
of  a  number  of  smaller  projectiles.  Furthermore,  the  hyper¬ 
velocity  projectiles  would  not  be  required  to  penetrate  the 
armored  vehicles  since  the  lethality  of  the  fragments  spalled 
from  the  rear  surface  of  the  armor  would  be  expected  to 
neutralize  the  soft  components  inside  the  vehicle.  For  example, 
projectiles  of  1  to  5  grams  with  velocities  in  the  neighborhood 
of  10  km/sec  would  cause  spallation  on  1-inch  armor  plate 
currently  in  use.  Of  course,  one  of  the  most  critical  items  in 
such  a  weapon  is  a  launcher  that  is  sufficiently  portable  and 
compact  that  it  may  be  delivered  ballistically  and  that  is  also 
adequate  to  launch  a  number  of  small  projectiles  to  the  high 
velocities  that  would  be  required.  Most  launching  systems 
capable  of  accelerating  integral  projectiles  to  10  km/sec  are 


too  massive  to  be  weaponized. 


However, 


an  explosively-driven 


1 


launcher  (References  1-7)  capable  of  launching  multiple  1-  to 
5-gram  projectiles  to  10  km/sec  showed  promise  of  weaponization. 

The  objectives  of  the  six-month  program  reported  here  were: 
(1)  to  demonstrate  the  ability  of  an  explosive  driver  to 
launch  a  projectile  to  9  to  10  km/sec  within  a  total  length  of 

1  meter  and,  (2)  to  determine  the  feasibility  of  launching 
multiple  projectiles  using  the  above  launcher. 

The  first  objective  was  accomplished  by  launching  approxi¬ 
mately  25  projectiles  to  greater  than  9  km/sec.  Although  the 
launching  of  intact  multiple  projectiles  was  not  achieved,  such 
a  goal  is  believed  achievable. 


2 
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SECTION  2 
BACKGROUND 


The  program  performed  was  a  feasibility  study  to  evaluate 
the  possibility  of  developing  an  explosive  driver  hypervelocity 
launcher  capable  of  launching  simultaneously  many  small  pro¬ 
jectiles  to  velocities  of  10  km/sec.  Any  feasib’e  launcher 
developed  must  be  capable  of  being  incorporated  into  a  weapon 
system.  To  design  a  .launcher  system  effectively,  its  ultimate 
usage  must  be  kept  in  mind:  generation  of  rear  surface  spall  on 
armored  vehicle  targets.  The  following  subsections  present 
background  information  concerning  the  constraints  placed  on  the 
launcher,  rear  surface  spall  phenomena,  and  the  explosive 
driver  concept. 


2.1  CONCEPTUAL  SYSTEM  CONSTRAINTS 


The  parameters  shown  in  Table  1  define  goals  for  a  con¬ 
ceptual  system  (Reference  8) .  Note  that  the  restrictions 
placed  upon  the  feasibility  study  parameters  are,  in  most 
cases,  significantly  relaxed  from  the  conceptual  system 
constraints. 
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TABLE  I 

J 

'} 

HYPERVELOCITY 

LAUNCHER  PERFORMANCE 

PARAMETERS 

i 

Item 

Conceptual 
System  Parameters 

Feasibility 

Program  Parameters 

Impeller  mass 

2-50  kgm 

Less  than  .10  kgm 

l 

Total  warhead  length 
(include  acceleration 
portion) 

Less  than  1  meter 

No  restriction 

Impeller  diameter 

Less  than  40  cm 

No  restriction 

Environment 

Ground  level  Atm. 

No  restriction 

1 

Projectile  flight 
distance 

10-100  meters 

No  restriction 

{ 

I 

No.  projectiles  per 
warhead 

20-1000 

1-10 

t 

:  jj 

Mass  of  individual 
projectile 

0.5-5  gm 

1  gm 

Density  of  individual 
projectile 

0.7-20  gm/cm3 

To  be  determined 

1 

j 

Projectile  aerodynamic 
drag  coefficient 

Less  than  0.2 

No  restriction 

> 

Projectile  ejection 
velocity 

6-10  km/sec 

9-10  km/sec 

, 

Impact  velocity 

5-8  km/sec 

No  restriction 

• 

Ejection  half  angle 
(with  respect  to 
warhead  axis) 

15°-40° 

15°-20° 

; 

: 

Prime  target?  steel 
armor 

1.5-6  cm  thick 

No  restriction 

I 

» 
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Explosive  drivers  have  been  used  in  the  past  to  launch 
single  projectiles  of  different  shapes  to  velocities  of  up  to 
12  km/sec.  However,  these  launcher  systems  were  considerably 
longer  than  the  1-meter  length  required  for  the  conceptual 
system  given  in  Table  1.  These  launcher  systems  were  neces¬ 
sarily  long  to  have  sufficient  launcher  gas  mass  without 

* 

correspondingly  high  reservoir  pressures.  Therefore,  once  a 
required  system  volume,  length,  and  diameter  are  set  and  the 
areal  density  of  the  projectile  is  given,  the  minimum  required 
reservoir  or  initial  projectile  base  pressure  is  fixed.  Any 
relaxation  of  the  warhead  length  for  a  feasibility  study  would 
represent  a  reduction  in  the  reservoir  pressure  and,  therefore, 
failure  to  establish  operational  feasibility  at  the  required 
system  pressure. 

The  launcher  system  to  be  developed  would  therefore  have 
the  following  restrictions:  length — 1  meter;  diameter — 
less  than  40  cm;  muzzle  velocity — 9  to  10  km/sec.  All  of  the 
tests  would  be  conducted  at  range  pressures  of  1  atmosphere  and 
in  air.  By  adhering  to  these  constraints,  the  feasibility  of 
developing  a  weaponizable  system  for  launching  multiple  hyper¬ 
velocity  projectiles  could  best  be  evaluated. 


The  final  gas  pressure  in  the  reservoir  is  proportional  to  the 
initial  gas  density  and  the  square  of  the  explosive  detonation 
velocity. 


2.2  MECHANICS  OF  REAR  SURFACE  SPALL  INDUCED  BY  PROJECTILE 

IMPACT 

The  impact  of  a  high  velocity  projectile  upon  a  slab  of 
material  results  in  a  compressive  stress  pulse  in  the  target 
that  propagates  away  from  the  point  of  impact.  This  stress 
pulse  contains  most  of  the  momentum  carried  by  the  projectile 
prior  to  the  impact.  When  this  compressive  pulse  encounters  a 
free  surface,  it  reflects  and  produces  tension  near  the  surface; 
if  the  cohesive  strength  of  the  material  is  exceeded,  fracture 
occurs,  and  some  or  all  of  the  momentum  of  the  stress  pulse  is 
carried  off  in  the  resulting  fragments  of  the  material.  This 
process  of  dynamic  fracture  by  reflection  of  a  compressive 
stress  pulse  at  a  free  boundary  is  termed  "spall." 

The  character  and  depth  of  the  spall  are  greatly  influenced 
by  the  amplitude,  shape,  and  spatial  width  of  the  compressive 
stress  pulse.  These  factors  are,  in  turn,  dependent  upon  the 
size,  velocity,  and  composition  of  the  projectile,  and  the 
material  properties  of  the  target. 

To  illustrate  how  the  nature  of  the  compressive  pulse 
influences  the  characteristics  of  spall,  let  us  consider  two 
simple  cases.  First  assume  a  rectangular  elastic  pulse,  as  in 


Figure  1,  with  amplitude  o  ,  moving  toward  a  free  surface  at 

x  .  For  this  discussion  stress  will  be  treated  as  f'ositive  in 
o 

compression.  The  material  contained  within  the  pulse  is  moving 
toward  xQ  with  a  velocity  vQ  =  oQ/pc,  where  p  and  c  are  the 
density  and  sound  speed  of  the  mate 'ial  (Figure  1).  Elsewhere, 
the  'laterial  is  at  rest.  When  the  loading  portion  of  the 
•'uise  mesL  .  '-ue  free  rear  surface,  it  reflects  and  becomes  an 
uv"=  o.»di-a  wave,  relieving  the  stress  in  the  material  ahead  of 
it  i'*  a  stress-free  state,  but  in  the  process  increasing  its 
velocity  by  a  factor  of  two  (Figure  2) .  As  the  right-going  and 
%f*;_yoing  unloading  waves  cross,  they  each  enter  unstressed 
nii  .'■rial.  Each  wave  reduces  the  magnitude  of  the  stress  in 
the  material  ahead  of  it  by  an  amount  (jq,  so  a  tensile  stress  of 
.is  produced  as  in  Figure  3. 

’  l.  che  magnitude  of  the  tensile  stress  at  this  point  is 
suifj. 'ient  to  produce  fracture,  spall  will  occur,  and  the 
mater: ul  to  the  right  of  the  failure  location  will  move  off 
with  velocity  somewhat  less  than  2vo,  leaving  the  material 
behind  at  rest,  and  with  both  portions  stress  free  (or  nearly  so). 


(a)  Stress  versus  distance  (b)  Particle  velocity  versus  distance 


Figure  1  Rectangular  stress  pulse  prior  to  interaction  with 

free  surface  at  x  . 
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la)  Stress  versus  distance  (b)  Particle  velocity  versus  distance 


Figure  2  Rectangular  stress  pulse  reflecting  from  a  free 
surface  prior  to  the  onset  of  tension. 


(a)  Stress  versus  distance  (b)  Particle  velocity  versus  distance 


Figure  3  Onset  of  tension  for  a  reflected  rectangular  stress 
pulse. 
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Note  that  the  initial  tension  (and  hence  failure)  develops  at  a 
distance  from  the  free  surface  equal  to  half  the  spatial  width 
of  the  incident  compressive  stress  pulse.  If  no  fracture 
occurs,  the  right-going  unloading  wave  reflects  from  the  free 
surface  and  becomes  a  left-going  loading  wave  that  increases  the 
negative  stress  (tension)  to  an  unstressed  level  (Figure  4). 

Next,  consider  a  triangular  compressive  stress  pulse,  as  in 
Figure  5.  The  development  of  the  tension  is  somewhat  easier  to 
visualize  in  this  case  by  construction  of  a  virtual  reflected 
pulse  that  becomes  real  as  it  enters  the  material.  The  stress 
at  each  instant  is  then  the  superposition  of  the  portions  of  the 
two  pulses  that  lie  within  the  material.  The  stress  state  at  a 
later  time— when  the  reflected  pulse  enters  the  material-is 
shown  an  Figure  fi.  Note  that  the  largest  tension  at.  each 
instant  results  from  the  superposition  of  the  peak  of  the  left¬ 
going  wave  with  successive  portions  of  the  right-going  wave. 
Since  the  right-going  wave  is  compressive,  it  always  acts  to 
reduce  the  tension.  Once  the  left-going  wave  passes  tht  last 


a 


(a)  Stress  versus  distance  (b)  Particle  velocity  versus  distance 


Figure  4 


Stress  pulse  after  complete  reflection  by  a  free 
surface  at  x  .  1 

o 


portion  of  the  right-going  wave,  no  further  increase  in  tension 
occurs.  Since  both  waves  are  moving  with  the  same  speed  (or  nearly 
30,  even  when  non-linear  materials  are  encountered)  in  the 
absence  of  spall,  the  peak  tension  is  again  developed  at  a 
distance  from  the  free  surface  equal  to  half  the  spatial  width 
of  the  incident  stress  pulse. 

Now  let  us  consider  how  the  spall  process  develops  when  the 
ocaV  stress  for  a  triangular  pulse  is  large  compared  to  the 
rpa* '  strength  of  the  material.  If,  for  instance,  a  spall  plane 
suddenly  formed  at  the  point  of  peak  tension  at  the  time 
illustrated  in  Figure  6,  a  new  free  surface  would  be  created 
there;  no  tension  would  remain  in  the  material,  but  the 
remaining  material  would  still  contain  a  right-going  compressive 
pulse  that  would  reflect  from  the  new  free  surface  and  hence 
develop  additional  tension  (Figure  7) .  The  process  of  the 
buildup  of  tension  and  spall  then  occurs  repeatedly  until  the 
leading  edge  of  the  reflected  pulse  passes  the  trailing  edge  of 
the  incident  pulse,  and  no  further  tensile  stress  increases 
occur . 


Real  pulses 


Virtual  pulses 

Unreflected  incident 
pulse  remaining  after 
spall 


Jew  virtual  reflected  stress  pulse 


Jew  free  surface 


Figure  7  Stress  pulse  (real  and  virtual)  for  remaining  material 
after  spall  has  occurred. 
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A  few  generalizations  can  be  made  from  these  simplified 
analyses.  The  maximum  depth  of  spall  will  not  exceed  half  of 
the  spatial  width  of  the  incident  stress  pulse.  For  a 
triangular  pulse  with  a  peak  that  greatly  exceeds  the  spall 
strength,  multiple  spalls  will  occur  and  will  probably  produce  a 
spray  of  very  fine  fragments.  A  more  rectangular  (flat-topped) 
stress  pulse  will  produce  fewer  spall  layers  and  will  lead  to 
comparatively  large  spall  fragments.  (For  an  introduction  to 
stress  wave  effects,  see  References  9  to  11.  More  advanced 
treatments  are  found  in  References  12  and  13.) 


The  impact  conditions  (i.e.,  projectile  size,  geometric 
configuration  and  velocity)  will  establish  the  initial  amplitude 
and  duration  (and  hence  spatial  width)  of  the  compressive 
stress  pulse.  The  material  properties  of  the  target  will 
affect  the  form  of  the  pulse  as  it  propagates  toward  the  free 
rear  surface,  especially  if  the  pulse  must  travel  a  large 
distance  compared  to  its  width.  For  many  materials  impacted  to 
high  stress  levels,  the  decrease  in  compressibility  with 
increasing  pressure  causes  the  stress  pulse  to  become  roughly 
triangular  after  propagation.  However,  the  pressure-induced 
phase  transition  that  occurs  in  iron  and  mild  carbon  steels  at 
around  150  kbar  leads  to  rarefaction  shock  upon  unloading, 
which,  in  turn,  causes  a  nearly  flat-topped  stress  pulse  to 
develop.  Consequently,  the  character  of  the  spall  fragments  in 
these  materials  may  be  quite  different  than  that  in  materials  in 
which  no  phase  changes  are  produced. 


Typical  armors  used  in  combat  vehicles  (iron  base 
armors)  undergo  this  phase  transition  around  150  kbar.  The 
unloading  shock  caused  by  this  phenomenon,  when  subjected  to 
projectile  impact,  leads  to  the  creation  of  rear  surface 
spall  fragments.  The  goal  of  the  program  being  conducted  is 
to  study  the  feasibility  of  a  launcher  capable  of  accelerating 
multiple  projectiles  to  hypervelocities  that,  when  impacting 
armored  vehicles,  would  utilize  this  rear  surface  spall 
mechanism  to  neutralize  "soft  components"  within  the  vehicle. 
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2.3  THEORY  AND  OPERATION  OF  THE  LINEAR  EXPLOSIVE  DRIVER 
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The  linear  explosive  driver  represents  a  technique  whereby 
a  substantial  portion  of  the  chemical  energy  of  an  explosive  is 
converted  in  a  controlled  manner  to  the  kinetic  and  internal 
energy  of  a  gas.  Basically/  the  energy  densities  in  the  gas  are 
produced  by  a  strong  shock  generated  by  the  progressive  collapse 
of  a  tube.  The  collapse  of  the  tube  is  such  that  it  may  be 
represented  as  a  piston  propagating  into  a  gas.  The  model  used 
to  describe  the  ideal  operation  of  the  driver  is  quite  similar 
to  that  used  to  explain  the  basic  discontinuous  motion  produced 
by  a  piston  in  one-dimensional  gasdynamics. 

The  operational  characteristics  of  the  linear  explosive 
driver  are  shown  in  Figure  8.  A  thin-walled  metal  tube  (the 
pressure  tube)  containing  the  driver  gas  is  surrounded  by  a 
chemical  explosive.  After  a  detonation  is  initiated  in  the 
explosive,  a  detonation  wave  propagates  in  the  explosive  along 
the  outside  of  the  metal  tube.  The  pressure  behind  the  detona¬ 
tion  wave  accelerates  the  tube  wall  in  toward  the  axis,  sealing 
the  tube  and  forming  a  conical-shaped  piston  (Figure  8b) .  The 
velocity  of  the  piston  is  equal  to  the  detonation  velocity  of 
the  explosive  (D) .  The  motion  of  this  piston  generates  a  strong 
shock  wave  in  the  stationary  column  of  the  driver  gas.  If  the 
gas  behaves  ideally  (i.e.,  the  ratio  of  the  specific  heats,  y, 
of  the  gas  is  constant),  then  the  velocity  of  the  shock  wave,  S, 
is  (y  +  1)  D/2.  The  position-time  histories  of  the  piston  and 
shock  wave  are  shown  in  Figure  8c  for  an  ideal  driver  gas 
(y  =  5/3) .  These  trajectories  are  presented  in  the  dimension¬ 
less  coordinates: 


x  =  J  and  t  =  jp 


i- 
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(dimensionless  time) 


Explosive 


Collapsed 
pressure  tube 


Detonation  wave 


Shocked  driver  gas 


Shock  wave 


(b) 

Time  t  =  t,,  after  initiation  of  the  explosive 


Position-time  history  of  shock  and  detonation  wave 
Figure  8  Ideal  operation  of  the  linear  oxplosio  drive: 


where  x  and  t  are  the  distance  and  time  after  the  shock  wave 
begins  to  move  ahead  of  the  detonation  wave,  d  is  the  internal 
diameter  of  the  pressure  tube ,  and  D  is  the  detonation  velocity 
of  the  explosive.  The  use  of  these  coordinates  facilitates  the 
comparison  of  drivers  having  pressure  tubes  of  different  dimen¬ 
sions  or  utilizing  explosives  with  different  detonation  veloci¬ 
ties.  It  should  be  noted  that  in  this  coordinate  system  all 
slopes  are  normalized  with  respect  to  the  detonation  velocity  of 
the  explosive.  For  example,  in  Figure  8c  the  trajectory  of  the 
detonation  wave  has  a  slope  of  unity,  while  that  cf  the  shock 
wave  has  a  slope  of  four^thirds: 


/  S  _  y  +  1  .  _  5  \ 

\D  "  2  -  '  *  "  3/ 


The  thermodynamic  state  of  the  shocked  gas  (subscripts  1) 
is  described  by  the  following  relations: 


Pressure : 


PX  -  »o  (’-H  °2  <» 


Temperature : 


_  V 


a  (^i) 


(2) 


Density: 


■  m 


(3) 


Sound  speed: 


■i  -  11 » 


(4) 


where  pQ,  PQ,  and  Tq  are  the  initial  density,  pressure,  and 
temperature  of  the  unshocked  gas.  These  relations  assume  that 


f 
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the  initial  pressure  in  the  unshocked  gas  is  quite  small  com¬ 
pared  to  the  corresponding  quantities  in  the  shocked  gas.  It  is 
noted  that,  for  a  given  ideal  gas,  the  magnitude  of  each  of 
these  properties,  except  density,  is  a  function  of  the  piston 
velocity  {detonation  velocity) . 

Using  the  ideal  theory  for  reflection  of  a  strong  shock 
from  a  rigid  wall,  the  gasdynamic  conditions  in  the  reservoir 
may  be  calculated.  The  pressure,  temperature,  and  sound  speed, 
for  example,  are  found  from  the  relations: 


The  performance  of  linear  explosive  drivers  has  been  demon¬ 
strated  over  a  wide  range  of  experimental  parameters.  Internal 
diameters  of  pressure  tubes  made  of  copper,  steel,  lead,  and 
aluminum  have  ranged  from  1/4  inch  to  16  inches;  high-explosive 
weight  has  ranged  from  27  grams  to  9200  pounds;  the  explosive- 
to-pressure  tube  mass  ratio  has  been  varied  from  0.5  to  10; 
driver  gases  have  included  helium,  air,  argon,  and  hydrogen;  the 
initial  pressure  of  the  driver  has  been  varied  from  15  to  2450 
psi;  and  the  detonation  velocities  of  the  explosives  used 
(liquid  and  solid)  have  ranged  between  5.5  km/sec  and  8.6 
km/sec.  While  some  of  these  cents  have  been  one  of  a  kind,  the 
majority  were  essential  to  comprehensive  experimental  studies  of 
a  particular  explosive  driver  c.esign. 
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Although  the  performance  of  many  of  the  drivers  adheres 
closely  to  the  ideal  predictions,  certain  departures  may  occur. 
Four  phenomena  have  been  incorporated  in  a  model  of  explosive 
driver  operation  to  account  for  observed  departures  from  the 
ideal  driver  performance  described  above.  They  are  (1)  radial 
expansion  of  the  pressure  tube  behind  the  shock  wave,  (2)  decom¬ 
position  or  predetonation  of  the  driver  explosive  during  the 
period  of  pressure  tube  expansion,  (3)  the  effect  of  boundary- 
layer  growth  behind  the  shock  wave,  and  (4)  formation  of  a 
metal,  gas,  or  metal-gas  jet  by  the  collapsing  pressure  tube. 
These  phenomena  are  inter-related  through  the  kinetics  produced 
by  the  imploding  pressure  tuba.  Their  interdependence  is  such 
that  changes  in  driver  behavior  resulting  from  certain  experi¬ 
mental  parameter  changes  cannot  be  attributed  solely  to  a 
particular  phenomenon.  However,  the  ability  of  the  model  to 
explain,  predict,  and  control  the  behavior  of  explosive 
drivers  justifies  the  categorization  of  these  phenomena.  A 
detailed  discussion  of  these  four  phenomena  is  presented  in 
References  1  through  4  and  14. 

In  a  typical  gun  design  the  explosive  driver  is  coupled  to 
the  barrel  of  the  gun  by  a  massive  steel  reservoir  section. 
Materials  used  to  form  the  reservoir  section  during  the 
course  of  launcher  development  included  lead,  steel,  concrete, 
and  explosive.  Guns  have  been  operated  in  both  the  chambraged 
and  unchambraged  mode.  When  the  strong  shock  generated  by  the 
explosive  driver  reaches  the  chambrage  plane  (or  projectile 
location  in  an  unchambraged  gun),  it  reflects  and  forms  a 
reservoir  of  very  high  enthalpy  gas.  The  reservoir  of  gas  is 
then  expanded  to  accelerate  the  projectile,  as  illustrated  in 
the  example  of  a  chambraged  gun  in  Figure  9. 
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projectile 


Figure  9  Operation  of  a  single-stage  explosively- driven 
launcher. 


For  launcher  applications  the  length  of  shocked  gas 
generated  by  the  explosive  driver  is  usually  limited  to  less 
than  10  tube  diameters  to  preclude  the  major  effects  of  boundary- 
layer  growth  or  pressure  tube  expansion.  The  jetting  of  the  col¬ 
lapsing  pressure  tube  that  forms  the  piston  is  usually  neg¬ 
ligible  in  launcher  applications  because  of  the  high  gas 
pressures. 


SECTICN  3 

EXPERIMENTAL  PROGRAM 


The  program,  as  conceived,  was  essentially  an  experimental 
program  designed  to  provide  sufficient  information  to  evaluate 
the  feasibility  of  the  successful  development  of  a  multiple  pro¬ 
jectile  hypervelocity  launcher  utilizing  explosive  drivers.  In 
this  role  the  program  was  structured  as  a  "cut-and-try"  evolu¬ 
tion  of  possible  launcher  designs.  Successful  features  of  launchers 
were  incorporated  into  subsequent  designs, and  features  deemed  to 
be  disadvantageous  were  replaced  with  more  probable  designs. 
Throughout  the  program  the  design  of  the  launcher  system  was 
aided  by  computational  ana?.ysii>.  This  computational  effort 
yielded  insight  into  the  design  of  the  very  high  pressure 
explosive  driver  system,  and,  after  the  successful  development 
of  this  driver,  aided  in  attempts  to  evolve  a  launch  cycle 
capable  of  accelerating  intact  projectiles  to  the  required 
10  km/sec. 

In  Lne  following  section  the  rationale  for  tse  original 
design  parameters  of  the  launcher  system  will  be  developed. 

The  results  of  the  testing  program  on  actual  launcher  hardware 
will  then  be  presented. 

3.1  INITIAL  DESIGN  OF  LAUNCHER  SYSTEM 

Past  attempts  to  launch  individual  projectiles  to 
hypervelocities  have  been  more  successful  with  projectiles 
constructed  of  an  alloy  of  magnesium  and  lithium,  or 
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Maglith.  This  metal  has  a  density  of  1.39  gm/cm  and  a  yield 
strength  of  about  20,000  PSI.  It  was  decided  in  the  initial 
conceptualization  of  the  launcher  that  Maglith  projectiles 
would  be  used  for  this  feasibility  study  also. 

A  1-gram  sphere  of  Maglith  has  a  diameter  of  1.11  cm. 

Allowing  for  an  additional  C.15  cm  thickness  of  Maglith  for 

a  sabot,  the  areal  density  of  the  launch  package  would  be 
2 

2.0  gm/cm  .  The  areal  density  was  thus  fixed  for  the  program. 
Even  if  projectile  materials  were  changed,  the  overall  areal 
density  of  the  projectile/sabot  package  would  have  to  remain 
the  same.  If  the  areal  density  were  to  increase,  a  longer 
barrel  would  be  needed  to  accelerate  the  assembly  to  a  given 
velocity  with  the  same  pressure  on  the  base  of  the  projectile. 
Figure  10  is  an  example  of  a  projectile/sabot  assembly  suitable 
for  launching  multiple  projectiles  from  one  barrel  while  main¬ 
taining  a  relatively  low  areal  density. 

Before  one  can  launch  a  projectile/sabot  assembly  such  as 

illustrated  in  Figure  10,  one  must  first  be  able  to  launch  a  flat 

2 

plate  of  areal  density  equal  to  2.0  gm/cm  .  To  reach  the  final 
velocity  for  which  a  launcher  system  is  designed,  the  projectile 
must  be  intact  throughout  the  launch  cycle.  If  the  projectile, 
or  a  projectile/sabot  assembly,  breaks  up  during  the  launch 
cycle,  then  the  velocity  would  be  less  than  that  achieved  by 
an  intact  projectile.  It  is  obvious  that  if  the  projectile/ 
sabot  assembly  were  to  break  up  in  the  barrel,  the  integrity 
and  velocity  of  the  individual  projectiles  would  be  questionable 
at  best. 

An  overall  length  of  about  1  meter  was  selected  as  the 
basic  requirement  for  the  driver  and  barrel  system.  Basic  physi¬ 
cal  principles  show  that  very  high  pressures  would  be  required 
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Figure  10  Possible  projectile-sabot  assembly  and  launch  sequence. 
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to  launch  projectiles  to  10  km/soc  in  such  short  lengths.  Since 
the  1  meter  length  included  both  the  explosive  driver  section 
and  the  barrel  length  required  to  reach  10  km/sec,  the  necessity 
of  high  chamber  pressure  is  even  more  clear. 

Experience  wich  the  explosive  driver  has  shown  that  the 
most  effective  method  of  creating  and  containing  very  high 
pressure  gases  is  through  the  use  of  a  linear  explosive  driver. 
Alternate  explosive  driver  techniques,  such  as  an  inverted 
driver  system,  were  considered  but  were  rejected  on  the  basis 
of  the  very  high  base  pressures  that  must  be  generated  and 
contained. 


For  design  purposes  it  was  assumed  that  the  launch  cycle 
could  be  approximated  by  the  launch  cycle  calculated  in  "The 
Theory  of  High  Speed  Guns,"  by  A.  Siegel  (Reference  15).  Past 
experience  with  explosive  driver  guns  has  shown  this  approxima¬ 
tion  to  be  relatively  accurate;  some  inherent  system  losses  are 
overcome  by  the  collapse  of  the  reservoir  from  external,  high 
explosive  tamping. 


In  the  design  calculations,  the  configuration  of  the 
launcher  is  assumed  to  be  that  of  Figure  11.  A  closed 
breech  was  assumed  which  contained  a  static  gas  of  mass  (G) , 
an  initial  pressure  of  (P^) ,  and  an  initial  sound  speed  of  (a2) . 
As  the  gas  expands,  it  accelerates  a  piston  of  mass  (M)  and 
area  (A) ,  The  assumption  made  in  the  gas  dynamics  cycle  is 
that  the  pressure  (P)  of  any  element  of  gas  during  the  expansion 
can  be  given  by 

P  =  P2  (V2/V)V 


?*!^$$H!?*V-W[y-!  ”»3S^  *  * ' y  **1? 


where 


Vj  is  the  initial  rolume  of  that  gas  element 

V  is  the  expanded  volume  of  that  gas  element 

Y  is  the  ratio  of  specific  heats  for  the  gas  and  is 
assumed  to  be  constant  {for  helium,  y  =  5/3). 


Figure  11  shows  the  curves  (from  Reference  15)  which  relate 
the  projectile  velocity  (Up)  and  the  required  barrel  length 
(Xp)  when  the  above  parameters  arc  known.  Perusal  of  these 
curves  indicates  that  high  projectile  velocities  are  facilitated 
by  having  as  high  a  sound  speed  as  possible  in  the  gas.  This 
sound  speed  is  proportional  to  the  detonation  velocity  of  the 
explosive  in  the  driver  {see  Equation  7,  Section  2.3).  For  the 
system  design,  therefore,  it  was  assumed  that  tne  explosive  used 
would  be  similar  to  an  HMX-based  plastic  bonded  explosive  (PBX) 
of  the  type  developed  and  used  by  the  Atomic  Energy  Commission. 

A  typical  detonation  velocity  of  these  explosives  is  8.5  km/sec 
with  a  density  of  1.87  gm/cc. 


In  the  geometry  of  Figure  9,  it  is  assumed  that  the  gas 
from  the  driver  stagnates  when  the  shock  impacts  and  reflects 
from  the  projectile/sabot  assembly.  Under  this  assumption,  the 
conditions  of  the  gas  can  be  calculated  to  a  close  approxima¬ 
tion.  After  trial  and  error  it  was  decided  that  a  25-kbar  driver 
should  be  used;  i.e.,  the  gas  behind  the  initial  shock  wave  is  at  a 
pressure  of  (P^)  of  25  kbar.  The  stagnated  pressure  (P2)  would 
be  (from  Equation  5) : 


P2  =  ePj^  =  150  kbar 


The  sound  speed  in  the  stagnated  helium  would  be  (from 
Equation  7) 

a2  =  9.8  km/sec 

2 

As  previously  stated  M/A  was  fixed  at  2.0  gm/cm  .  Reference  to 
Figure  11  shows  that  a  projectile  velocity  of  9  km/sec  (i.e., 
Up/a2  =  0.92)  can  be  achieved  with  G/M  ratio  of  1.  Then  from 
Figure  11,  the  barrel  length  was  determined  to  be 


The  initial  gas  density  required  in  the  pressure  tube  (pQ) 
of  a  25  kbar  driver  can  be  obtained  from  the  Equation  1: 

2P1 

p  =  - =•  =  0.026  gm/cc 

°  (y+1)D^ 

The  required  areal  density  of  the  gas  (Sg)  for  a  given  gas  mass 
to  projectile  mass  ratio  is 

6g  =  (§)  6p  =  (1.0) (2.0)  =  2.0  gm/cm2 

The  initial  length  of  driver  gas  (i^)  is  then 

JL  ==  ^  =  77  cm 
9  Po 

The  combined  length  of  barrel  and  driver  is  then  99  cm, 
which  is  just  within  the  specification  of  1  ueter  overall 
length.  The  length  could  have  been  shortened  by  using  a  higher 
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pressure  driver.  It  was  felt,  however,  that  the  launching  pres¬ 
sures  should  be  kept  as  low  as  possible  to  insure  the  integrity 

of  the  projectiles  during  the  launch  cycle.  Additionally,  a 

3 

pressure  of  2100  PSI  is  required  to  obtain  0.026  gm/cm  helium 
density.  This  pressure  is  almost  the  limit  that  can  be  obtained 
from  conventional  gas  bottles. 

A  word  is  in  order  concerning  the  temperature  of  the  helium 
gas  and  the  assumption  of  a  known  and  constant  y .  Figure  12 
shows  the  calculated  temperature  and  the  degree  of  ionization  at 
various  pressures  for  helium  gas  which  is  first  shocked  in  an 
explosive  driver  and  then  stagnated.  This  calculation  comes 
from  the  Saha  equation  (Reference  16).  Temperatures  are 
shown  assuming  both  nitromethane  and  PBX  to  be  the  driver 
explosive.  At  these  temperatures  ablation  and  radiation  effects 
upon  the  projectile  can  probably  be  neglected.  Many  projectiles 
of  plastic  and  metal  have  been  launched  at  Physics  International 
by  helium  at  this  temperature.  Flash  radiographs  of  the  pro¬ 
jectiles  have  shown  no  detectable  ablation  of  the  projectile  due 
to  the  launching  gas. 

Fro.;,  these  calculations,  the  basic  driver  and  launcher 
parameters  were  defined  and  shown  to  be  feasible.  An  explosive 
driver  that  generated  a  25~kbar  incident  shock  in  helium  gas  was 

required  to  launch  a  projectile  with  an  areal  density  of  2.0 

2 

gm/cm  .  To  achieve  these  pressures,  an  explosive  with  a  detona¬ 
tion  velocity  of  8.5  cm/psec  would  be  used,  with  the  helium  gas 
at  an  initial  pressure  of  2100  psi.  The  length  of  driver 
required  to  yield  a  G/M  =  1  would  be  77  cm,  and  the  barrel  for 
accelerating  the  projectile  would  be  22  cm  in  length. 
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Figure  12  SAHA  equation  for  helium. 
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With  these  basic  parameters  and  requirements  in  mind,  the 
experimental  program  was  begun.  Since  a  25  kbar  driver  had 
never  been  built  before,  this  was  a  totally  new  regime  in  which 
to  expex'iment.  Reference  was  made  to  previous  drivers  which  had 
achieved  15-kbar  incident  shocks,  and  from  the  results  obtained 
in  these  tests,  inferences  were  made  about  the  necessary 
features  of  a  25-kbar  driver. 

3.2  LAUNCHER  DEVELOPMENT  AND  TESTING 

Eight  experiments  i#ere  conducted  in  this  program  over  a 
period  of  six  months.  Three  of  these  experiments  were  attempts 
to  launch  multiple  projectiles.  The  other  shots  were  attempts 
to  launch  intact  flat-plate  projectiles  to  hypervelocities. 

These  flat  plate  projectiles  simulated  possible  projectile/sabot 
assemblies  that  would  launch  multiple  projectiles  from  a  single 
barrel.  Each  of  these  shots  will  be  discussed  in  the  following 
section.  Table  2  is  a  summary  of  the  shots  conducted  and  the 
results  of  the  experiment. 

The  first  test  conducted  in  the  program  was  designed  to  be 
a  checkout  of  the  explosive  driver  configuration  to  be  used  in 
the  later  phases  of  the  program.  An  aluminum  disk,  4.14  cm  in 
diameter  and  0.74  cm  thick,  was  used  to  simulate  a  seven-pro¬ 
jectile  sabot  assembly  to  be  utilized  in  later  experiments.  The 
thickness  of  the  disk  was  adjusted  to  yield  an  areal  density 
equivalent  to  the  anticipated  projectile/sabot  assembly. 

Since  working  with  a  liquid  explosive  is  simpler  than  hand 
packing  a  plastic  explosive,  nitromethane  was  chosen  for  the 
driver  explosive  in  the  first  shot.  The  detonation  velocity 
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Explosive  Driver  Characteristics 


Shot 

High  Explosive 

Observed 

Detonation 

Velocity 

(cra/Usec) 

Observed 

Shock 

Velocity 

(cm/usec) 

Resultant 

Incident 

Chock 

Pressure 

(kbar) 

Driver 

Gar. 

Blunderbuss 

-1 

Nitromethane 

0.73 

0.78 

11.8 

Helium 

Blunderbuss 

-2 

Composition 

C-4 

0.917 

1.055 

19.8 

Helium 

Blunderbuss 

-3 

Composition 

C-4 

0.526 

1.10 

21.5 

Helium 

Blunderbuss 

-4 

Composition 

C-4 

0.52 

1.01 

18.1 

Helium 

Blunderbuss 

-5 

Composition 

C-4 

0.921 

1.042 

19.3 

Helium 

Blunderbuss 

-6 

Nltrome thane 

0.675 

Data 

inconclusive 

Not 

retrievable 
Designed  to 
be  6.45 

Helium 

Blunderbuss 

-7 

Composition 

C-4 

0.916 

1.11 

21.9 

Helium 
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Projectile  Characteristics 


Target  Characteristics 


Dian 

(cm) 

4.14 


3.84 


3.94 


Projectile 

Thickness 

(cm) 

0.74 


0.74 


0.74 


Mass 

(gm) 

26.9 


23.1 


24.4 


Barrel 

Length 

(cm) 

20.0 


20.1 


20.0 


Number 

of 

Barrels 


Observed 

Launch 

Velocity 

(cm/usec) 

0.65 


Fragments 
0.848 
Main  body 
0.71 


Fastest 

Fragment 

0.91 

Large 

Fragments 

0.31 


Projectile 

Condition 

Down-Range 

Edge  of  projectile 
sheared  off. 
Projectile  warped  in 
flash  X-rays 

Broken  up  severely 
around  edges 
X-rays  do  not  clearly 
show  main  body  of 
projectile 

Completely  broken  up 
into  many  small 
fragments  and  some 
larger  fragments 


Target 

None 


1  inch  thick 
armor  plate 


1  inch  thick 
armor  plate 


2 

3 

1.44 

1.75 

3.25 

18.1 

2 

Fastest 

All  projectiles  broke 

1  inch  thick 

2 

1.44 

0.74 

3.25 

18.3 

3 

Fragment 

up 

armor  plate 

1.44 

0.257 

3.26 

18.3 

2 

0.6 

1 

3.73 

1.43 

21.9 

25.4 

1 

0.909 

Relatively  Intact 

1  inch  thick 

edges  appear  to  have 
been  broken 

armor  plate 

4 

1.27 

1.43 

2.54 

100 

4 

0.<79 

Intact 

Ik  inch  thick 

4 

1.27 

0.74 

2.53 

100 

4 

max 

0.390 

armor  plate 

4 

max 

1.27 

1.76 

2.54 

100 

4 

0.492 

max 

1 

3.66 

1.43 

21.1 

25.4  cm 

1 

0.953 

Badly  broken  up  into 

2  one-inch 

plus  12.7 

two  large  pieces  and 

thick  armor 

cm  slotted 

many  smaller 

plates 

barrel 

separated  by 

extension 

one  inch 

7 

1.09 

1.24 

1.62 

25.4  cm 

1 

0.849 

Totally  broken  up. 

1  inch  thick 

Sabot 

Sabot 

each 

plus  12.7 

Fragments  burned  up 

armor  plate 

diam 

thickness 

Total 

cm  unslotted 

in  atmosphere  as  they 

3.73 

1.43 

21.9 

barrel 

progressed  down  range 

extension 


Damage 

N/A 


Punched 
5  cm  dia  hole 
through  plat# 
with  severe 
rear  surface  spall 

Virtually 
no  damage 


virtually 

no  damage 


Severe  damage 
5  cm  dia  hole 
with  much 
rear  surface 
spall 

A1  6  Mg-Li-cratered 
plate  about 
1.3*  dia  »  0.25"  deep 
nylon  did  very 
little  damage 
no  spall  from 
any  impact 

Severe  damage 
to  first  plate, 
with  large  spall 
fragment  thrown 
off.  Nc  damage 
to  second  plate 

Virtually 
no  dwage 
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required  to  yield  a  25  kbar  driver  with  a  reasonable  (2100  psi) 
initial  gas  loading  pressure  is  0.85  cm/psec,  however,  nitro- 
methane  has  a  detonation  velocity  of  only  0.625  cm/psec.  To 
overcome  this  low  detonation  velocity,  the  phenomenon  of 
predetonation  of  nitromethane  was  to  be  used  to  increase  the 
detonation  velocity.  Reference  14  discusses  the  mechanism  of 
predetonation  in  nitromethane. 

As  shown,  in  Figure  13,  a  heavy  steel  tamper  was  placed 
around  the  nitromethane.  As  the  shock  wave  in  the  helium 
"breaks  out"  in  front  of  the  detonation  of  the  nitromethane,  it 
sends  a  shock  wave  through  the  undetonated  nitromethane.  When 
this  shock  reflects  off  the  massive  tamper,  very  high  pressures 
are  generated.  These  would  hopefully  predetonate  the  nitro¬ 
methane.  Past  experiments  conducted  at  PI  had  shown  such  an 
effect  with  drivers  designed  to  yield  4  kbar  shocked  gas, 
although  the  induction  time  for  this  effect  to  occur  was 
rather  long.  It  was  felt  that  since  the  experiment  would 
yield  at  least  a  15  kbar  shock  (based  on  the  0.625  cm/psec 
detonation  velocity  of  nitromethane) ,  the  induction  time 
required  for  the  predetonation  of  the  explosive  would  be 
considerably  reduced. 

Analysis  of  the  shot  data  shows  that  the  detonation 
velocity  of  the  nitromethane  did,  in  fact,  begin  to  accelerate; 
however,  adverse  effects  from  gas  loss  caused  by  incomplete 
closure  of  the  pressure  tube  prevented  the  explosive  pre- 
initiation  from  accelerating  to  the  detonation  velocity 
required.  Figure  14  is  a  plot  of  the  detonation  and  helium  gas 
shock  x-t  history.  Ion  pins  and  cap  pins  were  used  in  the 
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Figure  13  Configuration  of  the  first  experiment. 


Figure  14  Blunderbuss-1  x-t  plot 


diagnostics  to  observe  the  detonation  and  helium  shock 
trajectories.  Ion  pins  are  simple  coaxial  conductors  with  a 
solid  sheath  that  shorts  out  when  the  conductive  ionized 
detonation  wave  passes  over  them.  Cap  pins,  on  the  other  hand, 
are  basically  the  same  as  ion  pins  with  the  addition  of  a 
metal  cap  over  the  end  of  the  conductor  which  is  sensitive  to 
pressure  and  shorts  out  the  coaxial  conductor  with  the  applica¬ 
tion  of  pressure.  A  simple  resistor-capacitor  discharge  circuit 
operating  at  ’’CO  volts  assures  high  signal  levels  from  the  pin 
to  the  oscilloscope.  Both  ion  and  cap  pins  were  located  every 
5  cm  along  the  length  of  the  pressure  tube.  Cap  pins  are 
passed  through  the  tamper  and  explosive  and  are  in  direct 
physical  contact  with  the  pressure  tube.  Ion  pins  are  passed 
through  the  tamper  and  are  in  contact  with  the  explosive  (see 
Figure  15) .  Since  the  shock  wave  does  not  precede  the  detona¬ 
tion  wave  until  breakout,  both  pins  react  to  the  detonation  wave 
up  to  that  point.  The  shock'  velocity  achieved  in  this 
experiment  was  0.78  cm/psec,  which  yielded  a  12-kbar  driver. 

The  flat  disk  projectile  was  launched  relatively  intact  to  a 
velocity  of  0.62  cm/ysec.  This  velocity  is  close  to  the 
ideal  velocity  that  would  be  achieved  for  a  12  kbar  reservoir 
pressure,  based  on  Seigel's  work.  The  fact  that  the  projectile 
was  launched  almost  intact  (a  thin  rim  section  of  the  disk  was 
sheared  and  was  launched  to  slightly  higher  velocities)  under  a 
12-kbar  reservoir  pressure  was  encouraging  and  pointed  out  that 
the  stresses  generated  by  the  launch  cycle  were  not  excessive  in 
terms  of  launching  an  intact  but  severely  distorted  projectile. 
Figure  16  is  a  flash  radiograph  showing  the  projectile  in 
flight. 
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Blunderbuss-1  was  unsuccessful,  partially  because  the 
predetonation  phenomenon  was  not  as  pronounced  as  was  hoped,  and 
partially  because  the  collapse  of  the  explosive  driver  pressure 
tube  was  incomplete.  Design  review  of  the  driver  indicated  that 
the  pressure  tube- to-explosive  mass  ratio  was  on  the  borderline 
of  full  collapse,  and,  from  the  results  of  the  test,  clearly  on 
the  wrong  side  of  the  line.  To  aid  in  the  further  design  of 
the  explosive  driver,  numerous  hydrodynamic  computer  code 
calculations  were  conducted.  The  results  of  these  calculations 
will  be  discussed,  in  Section  4.  The  most  important  reason  for 
conducting  these  calculations  was  to  aid  in  the  prediction  of 
the  pressure  tube- to- explosive  mass  ratio  that  would  insure 
complete  collapse  of  the  driver  pressure  tube. 

Since  it  appeared  that  even  under  very  high  pressures,  the 
predetonation  of  the  nitromethane  explosive  would  not  increase 
the  detonation  velocity  sufficiently  to  yield  a  25-kbar 
driver,  it  was  decided  that  an  explosive  of  a  higher  detonation 
velocity  should  be  utilized.  Various  explosives  were  con¬ 
sidered,  and  Composition  C-4  was  selected  as  the  candidate 
explosive. 

The  second  shot  of  the  program  was  again  a  driver  checkout 
shoe.  Figure  17  is  a  drawing  of  the  shot  configuration.  The 
driver  was  designed  to  yield  a  25-kbar  incident  shock  in  the 
helium  gas.  The  internal  pressure  of  the  helium  prior  to  the 
shot  was  2100  psi,  which,  when  coupled  with  a  detonation 
velocity  of  0.85  cm/usec,  would  yield  a  25  kbar  shock  in  the 
gas.  As  in  the  first  shot,  the  projectile  to  be  launched  was 
an  aluminum  disk  0.74  cm  thick  by  3.84  cm  in  diameter. 
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in  the  first  shot  the  edge  of  the  projectile  sheared  off, 
leading  to  partial  breakup  of  the  disk-like  projectile,  even 
though  the  reservoir  pressures  obtained  were  significantly  less 
than  the  pressures  required  to  achieve  10  km/sec  muzzle 
velocity.  To  eliminate  this  edge-shearing  phenomenon  in  the 
second  shot,  the  disk-like  projectile  was  designed  as  a 
truncated  cone,  with  the  edges  tapered  at  an  8  degree  angle. 
This  taper  would  seat  against  a  similar  taper  in  the  barrel 
liner  and  would  hold  against  the  initial  pressure  of  2100  psi. 


Previous  experience  using  Composition  C-4  in  explosive 
driver  systems  has  shown  a  tendency  for  the  detonation 
velocity  of  the  C-4  to  increase.  This  effect  occurs  due  to 
precompression  of  the  explosive  by  the  shock  in  the  helium  gas 
that  out-runs  the  detonation  in  the  explosive.  Thus,  while  the 
standard  detonation  velocity  of  C-4  is  0.804  cm/psec,  it  was 
felt  that  the  increase  in  detonation  velocity  due  to  precom¬ 
pression  would  approach  the  0.85  cm/psec  required  to  achieve  a 
25  kbar  driver.  To  insure  that  the  pressure  tube  would 
collapse,  efficiencies  and  C/M  ratios  typical  of  previously 
tested  high  pressure  drivers  were  used  in  this  new  driver 
system.  Additionally,  the  configuration  of  the  driver  system 
was  checked  out  on  the  1-1/2  D  ELK  code  that  was  assembled  to 
assist  in  the  design  of  the  drivers  for  this  program.  The 
results  of  this  calculation  showed  that  the  driver  pressure  tube 
would  indeed  collapse  completely. 


Examination  of  the  results  of  this  shot  show  some  very 
interesting  phenomena.  Initially,  the  detonation  velocity  of 
the  C-4  explosive  was  the  expected  0.804  cm/psec.  However, 
after  breakout  of  the  helium  shock  ahead  of  the  detonation 


41 


- — 


B&ffl5§@S5£fflfifflSEgiB8!ffi«s»ssssmwiBam*a«sMmm»»!am^ 


front,  the  detonation  velocity  of  the  C-4  increased  to  0.914 
cm/psec  due  to  the  precompression  of  the  high  explosive.  Fig¬ 
ure  18  is  an  x-t  plot  of  the  driver  for  this  shot.  This  increase 
in  detonation  velocity  was  far  above  the  expected  increase  and 
added  to  the  efficiency  of  the  driver  system.  However,  gas 
losses  due  to  bubble  entrapment  at  the  closure  point  of  t.he  collap¬ 
sing  pressure  tube  led  to  a  final  shock  velocity  of  1.055  cm/psec, 
cr  an  incident  shock  pressure  of  19.8  kbar.  Figure  19  shows  the 
collapsed  pressure  tube,  partially  sectioned,  releasing  areas  of 
gas  bubble  entrapment  in  the  collapsed  tube.  The  bubble  entrap¬ 
ment  is  a  manifestation  of  the  growth  of  boundary  layer  gases 
in  the  pressure  tube. 

The  projectile  was  launched  to  a  velocity  of  approximately 
7  kin/sec;  however,  it  was  not  an  intact  projectile.  Examination 
of  the  flash  radiographs  showed  that  the  edge  of  the  projectile 
had  broken  off  into  significant  pieces  and  that  various  small 
fragments  of  the  broken  pieces  were  accelerated  to  a  velocity  of 
8.5  km/sec.  In  the  radiographs  the  projectile  appeared  to  be 
tumbling;  however,  a  target  plate  of  1-inch-thick  armor  plate 
showed  a  very  neat  1.5-inch-diameter  hole  punched  through  the 
center,  indicating  that  the  projectile  impacted  the  plate  normal 
to  its  line  of  flight  (Figure  20). 

From  the  results  of  this  shot  it  was  concluded  that  a 
workable  20-  to  25-kbar  driver  had  been  designed  and  that  this 
design  would  be  utilized  in  future  experiments.  Although  the 
ultimate  incident  shock  obtained  by  the  driver  tested  achieved 
a  shock  of  only  20  kbar,  the  effort  required  to  eliminate 
the  gas  bubble  entrapment  and  increase  the  pressure  to 
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Figure  20a  Front  view  of  1  inch  thick  armor  target  after  impact. 


Figure  20b  Rear  view  of  1  inch  thick  armor  target  after  impact. 

Note  large  spall  area  compared  to  penetration  hole. 


25  kbar  was  not  deemed  necessary.  If  required,  a  small 
extension  of  the  barrel  length  would  compensate  for  the  loss 
in  reservoir  pressure. 

Evidence  from  the  previous  shot  indicated  that  the  edge 
of  the  projectile  was  breaking  up  from  shear  stre  -  induced 
when  the  projectile  began  acceleration.  In  the  previous  shot 
the  projectile  was  held  in  place  by  a  taper  designed  to  fail 
only  in  shear.  It  was  felt  that  a  new  design  that  allowed  the 
projectile  to  feel  very  little  shear  forces  at  the  beginning  of 
the  acceleration  would  increase  the  probability  of  launching  an 
intact  projectile. 

The  third  shot  of  the  experimental  series  utilized  the 
same  explosive  driver  as  that  developed  for  the  second  shot, 
Blunderbuss-2.  The  design  of  the  projectile  was  modified  so 
that  a  thin,  cylindrical  section  of  the  rear  of  the  projectile 
would  fail  in  tension  rather  than  allowing  the  edge  of  the 
projectile,  through  its  support,  to  fail  in  shear.  Figure  21 
is  a  drawing  of  the  shot  configuration.  In  this  configuration 
the  edges  of  the  projectile  would  not  see  shear  forces  induced 
by  the  breaking  up  of  the  projectile  support.  It  was  hoped  that 
this  design  would  yield  an  intact  projectile  at  the  required 
velocity. 

The  explosive  driver  of  Blunderbuss- 3  performed  con¬ 
sistently  with  the  previous  driver  performance,  yielding  an 
incident  shock  of  21.5  kbar.  However,  the  aluminum  disk-like 
projectile  was  completely  broken  up.  Very  small  fragments  of 
the  projectile  were  accelerated  to  a  velocity  of  9.05  km/sec. 
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Immediately  following  Blunderbuss-3,  Blunderbuss-4  was 
fired  which  was  designed  to  launch  seven  projectiles  to 
10  km/sec.  This  shot  utilized  projectiles  of  aluminum, 
nylon,  and  steel.  The  configuration  of  these  projectiles  was 
the  same  as  that  in  the  previous  shot,  i.e.,  they  were 
designed  to  fail  in  tension  along  a  cylindrical  parting 
plane  behind  the  projectile  itself  (Figure  22) .  This  shot 
utilized  seven  individual  barrels  to  launch  the  projectiles, 
each  projectile  being  approximately  0.5  inch  in  diameter 
(Figure  23).  To  insure  equal  acceleration  characteristics,  the 
areal  density  of  all  seven  projectiles  was  maintained  at 
2.0  gm/cm^. 

Diagnostics  on  this  shot  showed  that  the  driver  again 
performed  as  expected.  Once  more  the  incident  pressure  jump 
across  the  shock  wave  in  the  helium  gas  was  approximately 
20  kbar.  Flash  radiographs  on  this  shot  failed  to  operate 
correctly,  so  that  no  usable  data  was  gathered  from  these 
diagnostics.  A  target  of  1-inch-thick  armor  plate  was  used  as 
a  witness  plate  for  the  shot,  however,  and  it  showed  that  all 
seven  projectiles  had  broken  up  before  impacting  the  plate. 
Velocity  achieved  by  the  projectile  fragments  was  a  maximum  of 
6.05  km/sec. 
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Following  the  unsuccessful  Blunderbuss-4,  an  in-house 
review  of  the  technical  status  of  the  program  was  conducted. 
Experimental  evidence  had  shown  that  the  explosive  driver 
portion  of  the  launcher  system  was  working  adequately  and 
would  satisfy  the  requirements  of  the  program.  However,  no  shot 
was  completely  successful  in  launching  an  intact  projectile  to 
the  required  velocities.  In  closely  examining  the  conditions  of 
projectile  launch,  a  theory  that  explains  the  projectile  breakup 
in  the  early  stages  of  launch  was  developed. 
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0.567  in.  dia.  nominal* 


LL 


*Fit  to  individual  barrel 


Material 

A 

B 

r 

Number  required 

Nylon 

0.689 

0.130 

0.491 

2 

Aluminum 

0.291 

0.032 

0.550 

3 

Steel 

0.101 

0.032 

0.546 

2 

All  dimensions  are  in  inches 
Figure  22  niunderbuss-4 — projectile  design. 
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When  the  shock  wave  in  the  helium  gas  generated  by  the 
explosive  driver  stagnates  against  the  rear  of  the  projectile, 
it  experiences  a  six-fold  increase  in  pressure  (see  Section  2.3, 
Equation  5).  Thus,  with  an  explosive  driver  that  yields  an 
incident  shock  of  20  kbar,  the  reflected  (or  stagnation)  pres¬ 
sure  would  be  on  the  order  of  220  kbar.  At  early  stages  of  the 
launch  cycle,  therefore,  the  projectile  is  subjected  to  a  com¬ 
pressive  stress  of  approximately  120  kbar.  This  is  well  above 
the  threshold  for  plastic  flow,  and  the  projectile  would  attempt 
to  expand  radially  to  relieve  the  axial  compressive  stresses. 

In  a  tight  fitting  barrel,  which  all  of  the  shots  up  to  this 
point  had  had,  the  edges  of  the  projectile  would  be  expected  to 
push  outward  against  the  barrel.  In  hydrodynamic  flow  condi¬ 
tions,  the  difference  between  radial  and  axial  stress  (i.e.,  the 
deviatoric  stress)  would  be  expected  to  be  on  the  order  of 
10  kbar  or  less  (for  aluminum).  Thus,  the  projectile  would  be 
subjected  to  120  kbar  of  axial  compressive  stress  attempting  to 
accelerate  the  projectile  and  approximately  110  kbar  of 
radial  stress  resulting  in  very  large  drag  forces  at  the 
periphery  of  the  projectile.  The  shear  stresses  that  result 
from  interaction  of  these  forces  i3  extreme,  resulting  in 
breakup  of  the  projectile.  To  confirm  or  modify  this  theory,  a 
two-dimensional  computer  analysis  of  the  launch  cycle  conditions 
was  begun.  The  results  of  this  analysis  are  discussed  in 
Section  4. 


The  fifth  shot  of  the  experimental  program,  Blunderbuss-5, 
was  designed  to  prevent  the  induced  drag  forces  and  subsequent 
shear  stresses  from  affecting  the  initial  launch  conditions  of 
the  projectile.  To  accomplish  this  the  projectile  diameter  was 
made  undersize  with  respect  to  the  bore  of  the  barrel,  and  a 


51 


v  -  '*<'7 


layer  (approximately  0.015  inch)  of  Dow-Corning  silicone  grease 
was  introduced  between  the  projectile  and  the  barrel.  This 
layer  of  silicone  grease  would  place  a  strengthless,  low- 
friction  material  between  the  projectile  and  the  barrel  which 
hopefully  would  prevent  the  large  induced  shear  forces  from 
occurring. 

Additional  changes  in  the  configuration  of  the  launcher 
were  as  follows.  Almost  all  previous  experience  with  explosive 
driver  accelerators  has  been  in  launching  projectiles  of  magne¬ 
sium/lithium  alloy.  Since  very  little  experimental  data  has  been 
developed  in  launching  aluminum  projectiles  with  explosive  driver 
technology,  it  was  decided  that  a  return  to  Maglith  was  in  order. 
Also,  the  thickness  of  the  barrel  was  increased  to  the  point  that 
there  would  be  very  little  chance  that  premature  barrel  failure 
would  result  in  the  lowering  of  the  launch  velocity  through 
abnormal  gas  losses.  Figure  15  is  a  drawing  of  the  configuration 
of  Blunderbuss-5.  Figure  24  is  a  Polaroid  photograph  of  the  shot 
just  prior  to  firing. 

This  fifth  shot  was  judged  to  be  quite  successful.  The  ex¬ 
plosive  driver  worked  as  expected  and  yielded  an  incident  shock 
of  19.3  kbar.  Projectile  launch  velocity  was  increased  to 
9.1  km/sec.  Figure  25  is  a  flash  radiograph  of  the  projectile 
in  flight  at  the  velocity  of  9.1  km/sec.  It  is  obvious  that  the 
projectile  is  somewhat  broken  up;  however,  the  thickness  of  the 
launched  projectile  is  equivalent  to  the  original  thickness  of 
the  projectile.  Since  the  velocity  achieved  by  the  pr-  ,cctiie 
is  quite  high,  it  is  unlikely  that  the  projectile  broke  up  during 
the  launch  cycle. 
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It  was  theorized  that  since  the  projectile  could  not  have 
broken  up  in  the  barrel  and  still  have  achieved  the  high 
velocity  at  launch,  breakup  must  have  occurred  after  the  pro¬ 
jectile  exited  the  barrel.  Since  the  pressure  on  the  rear  of 
the  projectile  is  not  completely  relieved  by  the  time  launch 
occurs,  compressive  stresses  on  the  projectile  could  be  severe 
enough  to  cause  projectile  breakup  through  the  sudden  release  of 
the  radial  constrainment  of  the  projectile.  Since  the  flash 
radiograph  of  Blunderbuss-5  shows  that  the  thickness  of  the  pro¬ 
jectile  corresponds  to  the  original  thickness  of  the  projectile, 
yet  the  diameter  of  the  projectile  is  reduced  and  non-symmetric, 
it  appears  that  radial  spall  was  the  probable  mechanism  of  pro¬ 
jectile  breakup  after  launch. 

In  a  further  attempt  to  launch  an  intact  projectile,  the 
sixth  shot  of  the  series  (Blunderbuss-7)  utilized  a  barrel 
extension  to  relieve  the  gas  pressure  behind  the  projectile, 
yet  maintain  radial  constraint  on  the  projectile  itself. 

Figure  26  is  a  drawing  of  the  barrel  extension  that  wa3  inset 
and  welded  to  the  basic  launcher  configuration  of  Blunderbuss-5. 
For  this  shot  the  silicone  grease  of  Blunderbuss-5  was  replaced 
by  a  Teflon  sleeve  around  the  Maglith  projectile.  Figure  27  is 
a  drawing  of  the  projectile  assembly. 

This  shot  attained  marginal  success.  Driver  performs  ic<* 
was  once  more  consistent  with  previous  tests.  The  addition  c 
the  barrel  extension  and  the  Teflon  sleeve  did  increase  the 
launch  velocity  of  the  projectile  to  9.6  km/sec.  However,  the 
original  disk-like  projectile  was  broken  into  three  smaller 
fragments,  and,  from  the  report  of  the  flash  radiographs 
(Figure  28) ,  the  sum  of  the  parts  did  not  equal  the  original 
whole. 
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Some  calculational  effort  was  expended  at  this  point  to 
characterize  the  launch  process  and  to  examine  the  collapse  of 
the  reservoir.  These  calculations  matched  with  observed  experi¬ 
mental  results.  Indications  from  these  computations  were  that 
the  projectile  is  still  under  pressures  of  about  35  kbar  at  the 
time  that  it  exits  the  barrel.  It  appears  that  the  high  pres¬ 
sures  on  the  rear  of  the  projectile  were  great  enough  to  cause 
the  projectile  to  extrude  into  the  slots  cut  into  the  barrel 
extension  to  vent  the  gases.  When  the  now  cross-shaped  pro¬ 
jectile  comes  to  the  end  of  the  barrel  extension  slots,  with 
part  of  the  projectile  extruded  into  the  slots,  the  resulting 
impact  forces  cause  the  destruction  of  the  projectile. 

One  other  experiment  was  performed  during  this  time  frame. 
An  attempt  was  made  to  launch  twelve  projectiles  of  three  dif¬ 
ferent  materials  to  10  km/sec  from  a  6  kbar  explosive  driver 
launcher  system.  Unfortunately,  the  explosive  driver  section  of 
the  launcher  failed  to  collapse  completely,  and  thus,  the  pres¬ 
sures  required  to  reach  10  km/sec  were  not  obtained.  The 
average  velocity  of  the  projectiles  was  on  the  order  of 
4.5  km/sec.  From  the  review  of  the  shot  data  it  appears  that 
the  charge-to-mass  ratio  of  the  high  explosive  to  the  driver 
pressure  tube  was  i.-tlow  the  minimum  required  to  fully  collapse 
the  tube  against  the  6  kbar  shock  in  the  helium.  (Note  that 
this  was  a  totally  different  driver  and  launcher  design  than  the 
now  standard  1  meter.  Composition  C-4  launcher.)  The  design  of 
the  shot  had  been  checked  by  the  use  of  the  1-1/2  dimension 
computer  code  and  had  been  predicted  to  close  the  pressure  tube. 
The  closure  predicted,  however,  was  marginal.  In  future  designs 
it  is  apparent  that  a  greater  margin  of  collapse  energy  must  be 
used. 
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Blunderbuss-8 ,  the  final  shot  of  the  experimental  program, 
was  an  attempt  to  launch  seven  individual  projectiles  contained 
in  a  single  sabot  utilizing  a  single  barrel  launcher.  This  shot 
was  designed  to  show  that  multiple  projectiles  could  be  launched 
intact  from  a  20  kbar  driver  system.  The  sabot,  as  shown  in 
Figure  29,  was  designed  to  take  advantage  of  the  possible  radial 
breakup  of  the  projectile  as  it  exits  the  barrel.  If  the 
mechanism  of  projectile  breakup  in  the  two  previous  shots  was 
radial  spall  caused  by  the  sudden  unloading  of  the  projectile  as 
it  left  the  barrel,  this  same  mechanism  would  cause  the  sabot 
to  break  up,  yet  it  might  leave  the  individual  projectiles 
intact.  Additionally,  since  the  breakup  was  expected  to  occur 
radially,  this  would  aid  ir.  the  dispersion  of  the  projectiles. 

Figure  30  is  a  drawing  of  the  shot  configuration.  A  barrel 
extension  was  added  to  increase  the  length  of  the  barrel  by 
10  cm,  and  the  projectile  was  surrounded  by  a  layer  of  silicone 
grease,  as  in  Blunderbuss-5.  There  were  no  slots  in  the  barrel 
extension  for  this  shot,  its  function  being  to  increase  the 
muzzle  velocity  of  the  projectile.  The  projectile/sabot 
assembly  was  similar  to  the  previous  shots  in  that  it  was  a 
disk-like  package  containing  seven  individual  projectiles. 

Figure  31  illustrates  the  sabot  and  projectile  configuration. 

This  experimental  shot  was  essentially  a  calculated  risk. 

It  was  felt  that  the  causes  leading  to  projectile  breakup  were 
known.  Only  one  shot  remained  in  which  to  launch  multiple  pro¬ 
jectiles  and  exhibit  total  proof  of  concept ,  and  it  was  decided 
that  a  multiple  projectile  launch  should  be  attempted.  The 
shot,  however,  did  not  meet  expectations.  The  explo¬ 
sive  driver  system  worked  as  expected;  however,  the  projectile/ 
sabot  package  was  totally  broken  up,  and  a  velocity  of  only 
8.4  km/sec  was  obtained.  Projectile  and  sabot  fragments  shown 
in  the  flash  radiographs  were  tenuous  at  best  and  are  not  shown 
in  this  report  due  to  the  limited  definition  obtainable  in  half¬ 
tone  reproductions. 
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SECTION  4 

COMPUTATIONAL  SUPPORT 


4 . 1  INTRODUCTION 


The  computational  effort  associated  with  the  program  was 
directed  towards  two  basic  objectives.  The  first  objective  was 
to  provide  data  to  assist  in  designing  and  analyzing  the  per¬ 
formance  of  the  explosive  driver  used  as  the  energy  source  for 
launching  the  projectiles.  This  objective  was  achieved  by  per¬ 
forming  a  serie;:  jf  one  and  one-half  dimensional  (one¬ 
dimensional  Moving  Wall  Flume,  see  Appendix  B  for  a  description 
of  the  flume  option)  calculations.  The  second  objective  of  the 
computational  effort  was  to  perform  a  calculation  of  the  inter¬ 
actions  at  the  driver/barrel,  barrel/projectile,  and  driver/ 
projectile  interfaces,  and  the  response  of  the  projectile  to 
the  induced  stresses.  A  two-dimensional  calculation  of  the 
early  stages  of  the  launch  cycle  was  performed  to  obtain  this 
information. 

4.2  MATERIAL  MODELS  AND  EQUATIONS  OF  STATE 

The  yield  models  aad  equations  of  state  used  to  describe 
the  helium  driver  gas,  the  C-4  explosive,  and  the  steel  of  the  pres¬ 
sure  tube,  tamper,  and  barrel  were  identical  for  all  the  calcula¬ 
tions.  The  two-dimensional  launch  cycle  calculation  also  con¬ 
tained  a  magnesium-lithium  (Maglith)  projectile.  The  equation- 
of-state  models  used  for  these  materials  in  the  computational 
analyses  are  summarized  in  Table  III. 
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4.3  EXPLOSIVE  DRIVER  CALCULATION 

In  general  the  principal  information  to  be  obtained  from  an 
explosive  driver  calculation  is  whether  or  not  the  liner  col¬ 
lapses,  the  velocity  of  the  point  of  collapse,  the  peak  pressure 
at  the  shock  front  in  the  driver  gas,  and  the  velocity  of  this 
shock  front.  This  information  may  be  computed  most  efficiently 
using  a  one-dimensional,  explicit,  finite-difference  continuum 
mechanics  code  with  a  Moving  Wall  Flume  (variable  area  duct) 
capability.  Until  recently.  Pi's  explosive  driver  calculations 
of  this  type  were  performed  using  a  one-dimensional  Lagrangian 
code.  Experience  proved  this  technique  useful  but  indicated  tnat 
a.  more  precise  solution  of  the  state  of  the  driver  gas  could  be 
achieved  if  the  gas  were  described  as  an  Eulerian  fluid 
retaining  the  Lagrangian  description  of  the  pressure  tube, 
explosive,  etc.  A  one-dimensional  Eulerian  code,  PISCES  IDE, 
was  available  (Appendix  A),  and  the  PISCES  IDE  and  PISCES  1DL 
(Appendix  B)  codes  were  dynamically  coupled  in  a  manner  similar 
to  that  employed  in  Physics  International's  two-dimensional, 
coupled  Eulerian-Lagrangian  code,  ELK.  The  code  resulting  from 
the  coupling  of  the  PISCES  IDE  and  1DL  codes  was  used  to  perform 
the  driver  calculations  for  this  program.  In  these  calcula¬ 
tions  the  driver  gas  is  treated  in  an  Eulerian  frame  of 
reference  while  the  pressure  tube,  explosive,  and  tamper  are 
treated  using  the  conservation  equations  in  their  Lagrangian 
form. 


Figure  32  shows  the  computational  grid  for  the  final 
Composition  C-4  explosive  driver  configuration  at  a  time  of 
10  usee  after  explosive  initiation.  Also  shown  are  vector 
velocities  of  the  tamper,  explosive,  pressure  tube,  and  helium 
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at  this  early  time.  Figure  33  shows  a  portion  of  the  grid  and 
vector  velocities  for  this  computation  at  a  time  of  80  usee.  As 
can  be  seen,  the  pressure  tube  is  collapsing  on  axis  as 
required  and  by  a  time  of  80  nsec.  The  helium  shock  has 
broken  out  and  is  leading  the  explosive  detonation  front.  The 
calculational  results  show  the  velocity  of  the  intersection 
point  of  the  axis  and  pressure  tube  to  be  0.75  cm/usec.  The 
velocity  of  the  helium  shock  is  1.20  cm/usec.  The  peak  pres¬ 
sure  of  the  shocked  helium  s  25  kbar. 

This  calculation  and  the  corresponding  experiments  verified 
the  design  of  the  25- kbar  explosive  driver. 

4.4  TWO-DIMENSIONAL  LAUNCH-CYCLE  CALCULATION 


A  critical  aspect  of  any  gun  designed  to  launch  projectiles 
at  very  high  velocities  is  tc  insure  that  the  dynamic  loading  of 
the  projectile  does  not  result  in  projectile  breakup.  A  two- 
dimensional  calculation  aided  by  complimentary  one-dimensional 
calculations  was  performed  to  examine  the  loading  of  the  pro¬ 
jectile  for  the  launch  system  designed  in  this  program. 

A  complete  two-dimensional  calculation,  including  the 
explosive  driver  startup  and  complete  collapse,  is  a  very  complex 
and  lengthy  computation  and  was  beyond  the  scope  of  this  pro¬ 
gram.  Therefore,  2  one-dimensional  Moving  Wall  Flume  calcula¬ 
tions  were  performed  to  provide  a  reasonable  approximation  for 
the  time-history  of  the  pressure  tube  collapse  as  well  as  the 
final  shocked  gas  volume  and  a  velocity  versus  time  function 
describing  the  helium  gas.  Both  of  these  calculations  were 
identical  to  the  previously  described  explosive  driver  calcula¬ 
tion  except  that  the  full  tamper  and  barrel  of  the  final 
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launcher  design  were  included.  The  initial  finite-dif ference 
grid  for  these  calculations  is  shown  in  Figure  34.  The  differ¬ 
ences  between  these  two  calculations  was  solely  in  the  treatment 
of  the  right  boundary  condition  for  the  helium  gas  (Eulerian 
fluid) .  In  one  calculation  the  right  boundary  was  a  simple  flow¬ 
through  condition,  i.e.,  zero  velocity  gradient  at  the  boundary. 
In  the  second  calculation,  the  right  boundary  was  fixed,  i.e., 
zero  velocity  at  the  boundary.  These  cases  bound  the  real 
problem  of  an  accelerating  projectile  in  the  barrel.  The 
results  of  the  two  calculations  were  used  to  obtain  approximate 
descriptions  of : 

•  The  time-history  of  the  pressure  tube  collapse 

•  The  time-dependent  volume  and  velocity  field 
associated  with  the  helium  just  prior  to  interaction 
with  the  projectile 

These  results  were  used  as  initial  and  boundary  conditions  for 
the  two-dimensional  calculation  of  the  response  of  the  pro¬ 
jectile  during  the  first  10  usee  of  the  launch  cycle. 

Figure  35  shows  the  initial  zoning  for  the  two-dimensional 
launch  cycle  calculation.  The  non-shaded  region  represents  the 
helium  gas.  This  region  is  bounded  radially  by  a  fixed  boundary 
and  on  the  left  by  a  moving  piston.  The  prescription  for  the 
moving  piston  was  based  on  the  results  of  the  2  one-dimensional 
calculations  described  above  and  models  the  collapse  of  the 
pressure  tube.  The  zones  for  the  helium  gas  had  initial  veloci¬ 
ties  approximating  the  conditions  of  the  gas  behind  the  helium 
shock  front  as  determined  from  the  one-dimensional  calculations. 
The  initial  volume  of  helium  gas  was  also  based  on  the  one- 
dimensional  results.  The  shock  front  in  the  helium  was  at  the 


71 


Initial  zoning  fo: 


r» ^rc^r^oi3t«T.  a*»s«ssm!WBmB9M!W95M8ra^lM»®«^^ — 


base  of  the  projectile  at  time  t  =  0.  r.'he  zones  representing 
the  magnesium-lithium  (Maglith)  projectile  are  those  shaded 
zones  to  the  right  of  the  helium.  The  Maglith  zones  were 
initially  at  rest,  simulating  the  state  just  prior  to  the  pro¬ 
jectile  interacting  with  the  driver  gars.  The  zones  immediately 
above  the  helium  and  projectile  represent  the  steel  barrel.  The 
equations  of  state  and  yield  models  used  in  the  two-dimensional 
calculation  were  identical  to  those  previously  described  in 
Table  III. 

Figure  36  shows  the  computational  grid  at  the  final  time  of 
10  usee.  As  can  be  seen,  by  this  time  the  projectile  has 
separated  from  the  barrel  and  is  beginning  to  deform.  Figure  37 
shows  the  state  of  stress  in  the  projectile  at  times  of  2,  4,  6, 

0,  and  10  psec.  Until  the  time  of  about  6  psec,  the  gradient  of 
stress  in  the  radial  direction  (away  from  axis  of  symmetry)  is 
very  small.  However,  beyond  a  time  of  about  6  p sec  the  reflec¬ 
tion  of  the  stress  wave  from  the  free  surface  and  interaction  of 
the  projectile  with  the  barrel  result  in  severe  radial  stress 
gradients  being  established  in  the  projectile.  Between  a  time 
of  8  and  10  psec,  the  projectile  and  barrel  separate. 

Due  to  programmatic  restrictions,  the  zoning  used  in  the 
two-dimensional  launch  cycle  calculation  was  coarse,  especially 
with  respect  to  the  steel  barrel  in  the  region  near  the  pro¬ 
jectile.  This  had  the  effect  that  a  single  zone  in  the  barrel 
dominated  the  response  of  the  projectile.  This,  in  turn,  results 
in  a  somewhat  imprecise  determination  of  the  stress  field  m  the 
projectile.  A  much  more  finely  zoned  computation  would  be  required 
to  determine  the  details  of  the  time-dependent  stress  field  in  the 
projectile.  Therefore,  the  results  of  the  launch  cycle 
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calculation  trust  be  interpreted  with  caution.  However,  on  the 
basis  of  the  calculation  performed  it  is  considered  likely  that 
the  lauding  conditions  that  were  modeled  would  have  resulted  in 
projectile  breakup.  This  tentative  conclusion  is  consistent 
with  experimental  results  and  the  subsequent  introduction  of  a 
gap  between  the  periphery  of  the  projectile  and  the  barrel. 


SECTION  5 

ANALYSIS  AND  CONCLUSIONS 


Significant  accomplishments  have  been  made  during  this 
six  month  program  towards  the  goal  of  developing  a  multiple  pro¬ 
jectile  hypervelocity  launcher.  Although  the  actual  launch  of 
intact  multiple  projectiles  from  a  single  barreled  gun  was  not 
accomplished,  the  experimental  results  and  the  computational 
analysis  completed  during  this  program  give  added  confidence 
in  the  ultimate  feasibility  of  utilizing  an  explosive  driver 
in  a  weapons  system  to  launch  multiple  projectiles. 

5.1  ANALYSIS 

5.1.1  25-kbnr  Explosive  Driver.  Section  i.l  described  the 
requirements  in  performance  of  an  explosive  driver  that  could  be 
used  to  launch  projectiles  to  10  km/sec  in  a  20  cm  barrel  length. 
From  the  experiments  conducted,  it  is  evident  that  the  explosive 
driver  developed  in  this  program  is  sufficient  to  accomplish  the 
task. 

Examination  cf  the  performance  characteristics  of  the  25- 
kbar  explosive  driver  yields  some  interesting  insight  into  the 
operation  o>.  very  high  pressure  explosive  drivers.  Referring 
back  to  Figure  18,  the  x-t  diagram  from  shot  Blunderbuss-2,  it 
can  be  seen  that  neither  the  detonation  velocity  nor  the  shock 
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velocity  are  constant.  Up  until  the  time  of  breakout,  that 
point  at  which  the  helium  shock  begins  to  lead  the  detonation  in 
the  explosive,  the  ion  pins  and  cap  pins  report  at  the  same 
time.  From  Figure  18,  it  can  be  seen  that  the  detonation  velocity 
of  the  04  explosive  is  0.8  cm/usec,  the  detonation  velocity 
observed  in  detonation  velocity  tests. 

After  breakout  of  the  helium  shock  in  front  of  the  detona¬ 
tion  wave  occurs,  there  is  a  marked  acceleration  of  the  detona¬ 
tion  velocity  of  the  04  to  a  velocity  of  0.914  cm/wsec.  This 
is  an  increase  of  14  percent  over  the  conventional  detonation 
velocity.  The  mechanism  whereby  this  increase  in  detonation 
velocity  occurs  is  a  precompression  of  the  04  by  the  shock  in 
the  helium.  As  the  shock  in  the  helium  moves  out  in  front  of 
the  detonation  wave,  it  transmits  a  shock  through  the  metal 
pressure  tube  into  the  explosive.  For  the  driver  designed  here, 
the  pressure  transmitted  to  the  explosive  would  be  on  the  order 
of  25  kbar.  This  precompression  of  the  explosive  increases  it3 
density,  and  the  detonation  velocity  is  increased.  While  this 
phenomenon  increajes  the  energy  density  of  the  explosive,  the 
volume  of  the  explosive  is  decreased  due  to  the  precompression, 
and,  as  would  be  expected,  total  energy  is  conserved. 

As  the  detonation  velocity  increases,  the  shock  velocity  in 
the  helium  gas  also  increases  to  a  relatively  stable  velocity  of 
1.147  cm/usec.  This  shock  velocity  corresponds  to  an  incident 
shock  pressure  of  25  kbar.  However,  as  the  detonation  proceeds 
down  the  explosive  driver,  losses  begin  to  occur  that  degrade 
the  performance  of  the  driver.  The  most  significant  losses 
arise  from  the  entrapment  of  bubbles  of  high  pressure  gases  at 
the  collapse  point  of  the  pressure  tube  and  the  radial  expansion 
of  the  pressure  tube  behind  the  shock  front  due  to  the  extremely 
high  pressures.  This  loss  of  gas  behind  the  shock  front  sends 
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rare factions  towards  the  shock,  and  results  in  a  net  lowering 
of  the  incident  shock  velocity  and  a  corresponding  lowering  of 
the  ir.cidc.it  pressure. 

Again  from  Figure  18,  the  final  shock  velocity  is  about 
1.055  cm/usec.  This  equates  to  a  shock  pressure  of  19.8  kbar. 

A  reduction  in  shock  pressure  from  25  to  19.8  kbar  has  taken 
place  due  to  driver  losses. 

The  first  task  required  in  this  program  was  the  development 
of  an  explosive  driver  capable  of  producing  pressures  on  the 
order  of  25  kbar  in  order  to  launch  disk-like  projectiles  to 
velocities  of  10  km/sec.  An  explosive  driver  was  developed  that 
was  designed  to  yield  a  25-kbar  incident  shock  in  helium,  and 
which  actually  yields  an  incident  shock  on  the  order  of  20  kbar. 
Experimental  evidence  in  launching  projectiles  obtained  during 
this  program  indicates  that  the  driver  developed  is  sufficient 
for  the  program  requirements. 

5.1.2  Hypervelocity  Launcher  System.  This  program  was 
designed  to  study  the  basic  feasibility  of  a  method  that  would 
be  able  to  launch  multiple  projectiles  to  velocities  of 
>  9  km/sec.  As  such,  the  goals  of  the  program  were  not  to  design 
operational  or  prototype  systems  for  field  use,  but  rather  to 
examine  the  basic  explosive  driver  launcher  techniques  to 
determine  the  feasibility  of  developing  a  launcher  within  the 
parameters  of  the  conceptual  system  requirements. 

The  approach  to  the  problem  was  as  follows:  Two  capabilities 
must  be  developed  to  be  able  to  launch  multiple  projectiles  to 
hypervelocities  within  an  overall  length  approximating  one  meter. 

(1)  An  explosive  driver  must  be  developed  that  will  yield  reservoir 


2 

conditions  capable  of  launching  a  projectile  package  of  2  gm/cm 
areal  density,  ard  (2)  a  projectile-sabot  assembly  must  be  de¬ 
signed  to  launch  multiple  individual  projectiles  from  a  single 

2 

barrel.  In  order  to  maintain  a  2  gm/cm  areal  density,  this 
projectile-sabot  assembly  will  be  in  the  form  of  a  flat  disk. 
Therefore,  to  prove  feasibility  of  the  concept,  ;  method  of 
launching  this  disk-like  projectile  to  hyperveio.  ties  must 
first  be  developed. 

This  initial  feasibility  study  concentrated  on  the  ability 
to  launch  intact  disk-like  plates.  For  this  reason,  all  shots 
(with  the  exception  of  three)  dealt  with  the  development  of  a 
launch  cycle  capable  of  accelerating  a  flat  disk-like  projectile 
to  velocities  of  9  km/sec  (30,000  ft/sec). 

In  the  first  design  (Figure  13) ,  the  edges  of  the  plate 
(projectile)  broke  up  rather  early  in  the  acceleration  process. 
Several  other  methods  of  support  wore  tried  (Figures  17  and  21) 
without  success.  Intuition  (later  confirmed  by  a  two-dimensional 
calculation)  indicated  that  the  interface  between  the  steel 
barrel  and  the  outer  periphery  of  the  low-density  metal  plate 
was  the  source  of  the  problem.  At  least  two  problems  stem  from 
this  interface.  When  stresses  of  the  order  of  100  kbar  are 
applied  to  the  two  metals,  the  stress  across  the  interface  must 
also  be  of  this  magnitude,  '’’he  shear  stress  as  the  plate  tries 
to  slide  down  the  barrel  will  be  in  the  tens  of  kilobars.  The 
edge  drag  caused  by  these  stresses  can  result  in  the  edges 
lagging  the  rest  of  the  plate.  This  would  obviously  lead  to 
breakage  of  the  edges.  The  second  problem  comes  from  the  severe 
impedance  mismatch  between  the  two  metals.  This  mismatch  results 
ultimately  in  a  stress  gradient  across  the  interface  which 
causes  an  inward  radial  motion  of  the  periphery  of  the  plate. 
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With  axial  shock  waves  of  approximately  100  kbar,  the  radial 
motions  caused  by  such  phenomena  can  be  large  enough  to  cause 
break-up  of  the  plate. 

It  was  for  the  above  reasons  that  the  design  of  the  plate 
was  changed  to  that  shown  in  Figure  15.  The  space  between  the 
edge  of  the  plate  and  the  steel  barrel  was  filled  on  one  shot 
with  silicone  grease  and  in  another  with  Teflon.  It  was  hoped 
that  this  geometry  would  improve  the  impedance  mismatch  problems 
and  eliminate  the  edge  drag.  The  change  was  quite  effective— 
the  plate  must  have  held  together  until  the  plate  emerged  from 
the  barrel  in  order  to  achieve  the  observed  velocity.  Our  analysis 
indicates,  however,  that  at  this  time  the  plate  was  still  subjected 
to  approximately  35  kbar  of  axial  pressure.  The  edge  rarefactions 
resulting  from  the  sudden  release  of  such  a  high  pressure  would 
also  cause  a  plate  to  break-up.  In  the  final  shot,  PI  took  the 
chance  that  these  rarefactions  would  lead  to  the  expansion  and 
breakup  of  the  sabot  but  not  of  the  1  gram  particles  embedded 
in  the  sabot.  Unfortunately, both  the  sabot  and  the  particles 
broke  up. 


5.2  CONCLUSIONS 

The  objectives  of  this  six-month  program  were:  (1)  to  demon¬ 
strate  the  ability  of  an  explosive  driver  to  launch  a  projectile 
to  9  to  10  km/sec  within  a  total  length  of  1  meter  and,  (2)  to 
determine  the  feasibility  of  launching  multiple  projectiles  using 
the  above  launcher.  During  this  program  the  following  specific 
tasks  were  accomplished: 

•  An  explosive  driver  capable  of  producing  an  incident  shock 
in  helium  gas  of  20  kbar  has  been  produced  and  has  been 
shown  to  work  reliably  and  repeatedly. 


«  A  series  of  one-dimensional  and  one  and  one-half  dimen¬ 
sional  computer  calculations  was  performed  to  understand 
driver  performance  and  to  aid  in  the  two-dimensional  cal¬ 
culation  of  the  initial  stages  of  projectile  launch. 

•  Experimental  shots  were  conducted  with  a  1-meter-long 
launcher  that  were  successful  in  launching  25-gram  pro¬ 
jectiles  to  a  velocity  of  9.6  km/sec.  The  projectiles 
were  only  partially  intact  at  distances  of  2  feet  down 
range.  It  is  felt  that  the  cause  of  projectile  breakup 
is  interaction  between  the  projectile  and  the  barrel 
during  the  launch  cycle. 

•  A  two-dimensional  Lagrangian  computer  calculation  was 
performed  to  analyze  the  initial  launch  cycle  characteristics 
of  the  hypervelocity  launcher.  This  calculation  showed 

that  the  most  likely  cause  of  projectile  breakup  is  in¬ 
duced  inward  radial  velocities  in  the  projectile  due  to 
the  passage  of  the  shock  wave  in  the  steel  barrel  over 
the  projectile. 


With  the  results  of  the  computer  calculation  and  a  careful 
review  of  the  test  results,  it  is  felt  that  the  cause  of  projectile 
breakup  can  be  isolated  and  corrected.  If  this  deficiency  is 
corrected,  then  it  appears  feasible  that  a  launcher  can  be  pro¬ 
duced  using  the  techniques  developed  j.n  this  program  to  launch 
intact  multiple  projectiles. 
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SECTION  1 
INTKOIHICTION 


What  follows  is  a  description  of  a  numerical  method  for  the  solution 
of  the  time-dependent  equations  of  continuum  mechanics  expressed  in  an 
Eulerian  framework.  The  method  is  general  in  that  a  wide  range  of  mate¬ 
rial  models  may  be  incorporated  into  it.  The  differential  equations 
solved  are  the  Eulerian  counterpart  of  the  Lagrangian  equations  formulated 
by  Wilkins  (Reference  1).  The  difference  method  used  is  essentially  that 
made  popular  by  MacCormack  (Reference  2).  It  is  shown  to  be  of  second- 
order  accuracy  in  both  space  and  tisie  and  stability  criteria  arc  derived. 
The  results  of  several  sample  problems  using  this  method  are  presented 
and  compared  to  solutions  obtained  using  the  PISCES  1DL  computer  code. 
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An  equation  of  state  must  be  provided  to  determine  the  pressure,  p, 
as  a  function  of  the  relative  volume,  V=  p^/p,  and  the  Internal  energy, 

I-  -  PQe.  For  example,  a  typical  equation  of  state  might  be 

p  =  ax  (T|-  1)  +  a2(t)-  if  +  a3 <T|-  if  +  bt)E 

where  ax,a.,,  aa  and  b  are  constants  and  ri  *  1/V  *  p/pQ  . 

For  flume  symmetry  (d-4)  the  area  A(r,t)  must  be  supplied  (along 
with  initial  and  boundary  conditions)  to  complete  the  specification  of 
the  problem.  This  area  function  may  be  either  pre-specified  or  it  may 
be  determined  dynamically  during  the  course  of  the  calculation. 

For  the  other  geometries  (d =» 1,2,3)  we  need  additional  relations  to 
determine  Z^Cd*  1,2,3)  and  Eg(d=2,3)  as  these  are  not,  in  general,  deter¬ 
mined  completely  by  the  pressure,  p.  These  relations  are 

Er  =  -p  +  Sj 

ZQ  =  -P+sa 

where  the  stress  deviators  sx  and  s2  are  determined  from  the  following 
relations. 

The  first  stress  deviator  equation: 

fif(psi)  +  +  '^"r^  (8i+fM)  *  0  (d*  l*2»3) 

The  second  stress  deviator  equation: 

f^f»r)  +  ^<W«a)  +  §Mp|*  -  2^dH+  (d-l)^(s2  +  |u)  -0  (d-2,3) 


9a  15  ~  J  si  W*  1) 
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The  third  stress  dfViatoi:  equation: 


s.i“*£si  +  sJ  (d- 1,2,3) 
where  m  Is  the  sheer  modulus. 

The  stress  deviators  as  determined  above  are  then  sKXlified  by  a  yield 
condition.  For  example, 

(•?  +  se  +  4)  -  f  * 0 

where  Y0  is  tiie  material  strength. 

Next,  we  rewrite  the  above  equations  in  vector  form,  introducing  the 
variables 

P 

m  -  pu 

c  =  p(e  t  uz/2) 

S  ptj 
T  -  psa 

The  equations  assume  the  compact  form 
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SECTION  3 

THE  DIFFERENCE  EQUATIONS 


The  two-step  difference  method  used  to  approximate  the  solution  to 
Eq.  (1)  is  defined  by 

®;  *  *,■  -  £(?;,  -  »,-)  -  <^>vi  *p;,  •  w  -  * «;  <*•> 

«;**  -  £(?;-  ?;;,)  *  (=*-=/.,)  *  ^  v°j”]  <2b> 


where  the  subscripts  refer  to  a  spatial  mesh  of  points  with  spacing 
.'at  and  the  superscripts  refer  to  the  times 


‘n  ■  I 


where  At1  is  the  time  increment  that  the  solution  is  advanced  during 

cycle  i.  Also,  F.n  0.°,  and  Rn  equal  F(U,n),  G(U.n),  and  R(U,n), 
J  >  J  J  J  J  J 

respectively. 


Dj  is  the  artificial  viscosity  operator  given  by  the  expression 
„  M  i  n  r.i/^n  r;n\  i  n  n  i/^n  r?n  \ 

°J  -Cq/bT[,U>1'Uj,'V  J'  ‘  ,WJ’UJ-»H°J  “*WJ 

where  u"  is  the  material  velocity  at  position  Xj  and  at  time  tn  and 

C  is  a  constant  of  order  unity,  A  value  of  C  different  from  zero  is 
<1  7  q 

not  needed  for  .stability  in  calculations  for  which  the  linear  stability 
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analysis  applies,  but  may  be  needed  to  preserve  stability  when  computing  a 
strong  shock,  motion  near  the  axis  of  cylindrical  or  spherical  symmetries, 
or  other  violent  or  nonlinear  phenomena.  For  these  calculations,  a  value 
of  Cq  of  2.0  should  suffice  to  preserve  numerical  stability. 

The  finite-difference  method  of  Eq.  (2)  first  obtains  the  approximate 
value  U.  at  each  point  by  Eq.  (2a)  which  uses  forward  differences  to  ap¬ 
proximate  the  spatial  derivatives.  The  approximate  solution  is  then  used 
to  calculate  F^  ,  Gj  ,  and  which  are  used  with  backward  differences  in 
Eq.  (2b)  to  obtain  the  new  value  of  the  solution  uj1*1. 

The  specification  of  the  solution  is  complete  if,  in  addition  to  the 
prescription  (2)  for  advancing  the  solution  in  time,  we  provide  initial 
conditions  Uj°;  j  «  0,  1,  •••,  j^,  and  boundary  conditions  U^\  U  n  ; 

n  =  0  j  ...  n  Jmax 

’  ’  ’  max* 

The  accuracy  and  stability  of  the  difference  method  (2)  is  analyzed 
the  Appendix.  It  is  shown  there  that  the  method  la  of  second  order  accu¬ 
racy  and  that  the  linearized  stability  criterion  is  given  by  the  well 
known  Courant  Friedricks  Lewy  (CFL)  condition  ||(|u|  +  c)  s  1  where 
u  is  the  local  fluid  velocity  and  c  is  the  local  sonic  velocity.  The 
linearized  stability  criterion  is  a  necessary  but  not  sufficient  condition 
to  insu;e  stability  of  the  actual  difference  equations  which  are  nonlinear. 
Thus,  to  apply  the  CFL  criterion,  we  introduce  a  time  step  "safey  factor” 
(ts  (where  0  <  f{j  <  l)  and  determine  the  stable  time  step  from  the  con¬ 
dition 


Lt  r.  f 


Ax 


ts 


U  +c 


(3) 


The  CFL  stability  criterion  is  obtained  by  ignoring  the  effect  of  the 
vector  R.  It  is  easily  seen  that  the  Inclusion  of  R  does  not  significantly 
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affect  the  linearized  stability  criterion  quoted  above  as  the  change  in 
the  eigenvalues  01  the  associated  amplification  matrix  is  only  the  order 
of  At.  Also,  R  is  included  in  the  two-step  difference  scheme  in  such 
a  manner  that  the  method  retains  its  second-order  accuracy.  Although  the 
method  is  strictly  accurate  to  second  order,  the  vector  R  may  be  of  such 
a  magnitude  that  the  third-order  inaccuracy  in  R  dominates  the  other  error 
terms,  resulting  in  an  effective  loss  of  accuracy  in  the  spatial  deriva¬ 
tive  terms. 

To  insure  that  this  is  not  the  case,  ire  subject  the  tlsK  step  At 
to  the  additional  constraint 


At  * 


Ri 


(4) 


where  and  R^  a*e  the  first  components  of  the  vectors  V  and  R  and 

(fae*  is  *n  "accura -y  factor"  which  satisfies  0  £  f  <  1.  A  recom¬ 
mended  value  for  f  is  0.2. 

3C 

The  actual  time  step  At  used  in  each  cycle  of  the  calculation  is 
the  maximum  value  of  At  which  satisfies  the  two  criteria  of  Eqs.  (3) 
and  (4)  at  each  location  in  the  mesh. 
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SECTION  4  | 

SAMPLE  CALCULATIONS  ) 

t 

| 

The  results  of  several  staple  calculations  are  presented  in  graphical  j 

fora  in  this  section  for  the  purpose  of  illustrating  the  coaputational 
capability  of  the  numerical  aethod  presented.  Each  problea  was  run  both 
with  the  PISCES  ICE  code  and. with  the  PISCES  1DL  code  to  facilitate 
comparisons  between  the  accuracy  and  the  efficiency  of  the  two  codes  in 
solving  a. wide  class  of  one- diaans tonal  continuum  nechanlcs  problems. 

In  each  case,  the  two  code:,  used  the  same  number  of  xones  (or  aesh  points) 
and  were  run  for  the  saaa  amount  of  problea  time.  The  number  of  compu¬ 
tational  cycles  and  the  computer  tine  needed  to  run  differed  between  the 
two  codes  because  of  different  tiae-atep  stability  criteria  for  Lagranglan  ! 

and  Eulevian  codes.  Timing  teats  have  shown  the  zone  cycle  time  (the 
time  cv  compute  one  zone  for  one  cycle)  to  be  virtually  identical  in  | 

PISCES  1DL  and  PISCES  IDE.  1 

i 

\ 

Test  Problea  1;  Elastic-Plastic  Plane  Plate  Impact 

i 

At  time  t«0,  a  steel  plate  of  thickness  25  *x  (100  zones)  moving  at 
a  velocity  of  .004  ca/Msec  impacted  a  stationary  steel  plate  of  thickness 
25  cm.  The  steel  wes  modeled  with  p*7.85,  K»  1.69,  C«.S22,  and  a 
von  Mises  yield  strength  of  Yo*.003. 

The  results  of  the  Eulerian  and  the  Lagranglan  calculation  are  sum¬ 
marized  in  'figures  la  and  lb,  where  a  series  of  stress  profiles  at  dif¬ 
ferent  tines  are  superimposed.  No  artificial  viscosity  was  used  in  either 
calculation.  PISCES  IDE  required  102  cycles  to  run  this  problea  to  35  Msec 
while  PISCES  1DL  required  129  cycles. 


Praciiiif  w  Mail  />. » 


TCAH  73-14 


Test  Problem  2:  Flume  Test  A 

At  tine  t  *  0,  there  exists  s  gss  with  p»  .7:01,  y*5/3  in  a  straight 
flume  of  length  40  cn  <100  zones).  A  center  section  of  4  cn  initially 
has  an  energy  density  of  10  (hot  gas)  while  the  remainder  of  the  ga* 
has  an  initial  energy  density  of  1.5.Xl0“*  (cold  gas). 

Figures  2a  and  2b  show,  the  pressure  distributions  in  the  pipe  at 
three  successive  times,  the  IDL  problem  used' a  quadratic  artificial  vis¬ 
cosity  coefficient  of  2.0  and  a  linear  artificial  viscosity  coefficient 
of  0.1.  The  IDE  problem  used  CQ»2.0.  PISCES  IDE  required  115  cycles 
to  run  this  problem  to  250  yeec,  while  PISCES  1DL  took  354  cycles. 

Test  Problem  3;  Flume  Test  B 

This  problem  is  identicel  to  Test  Problem  2  except  for  the  shape  of 
tie  flume.  In  Test  Problem  2,  the  flume  was  straight.  In  this  problem, 
the  flume  has  a  straight  section  in  the  center  section  of  4  cm  and  in 
diverges  outward  from  the  center  each  way  with  a  slope  of  1/8. 

The  results  of  this  problem  are  Illustrated  in  Figure  3.  PISCES  IDE 
required  116  cycles  and  PISCES  1DL  requires  330  cycles  to  run  this  problem 
to  250  psec. 


Test  Problem  4:  Cylindrical  Implosion — Paynes  Problem 


This  is  s  problem  whose  solution  by  other  methods  appears  in  the 
literature  (References  5  and  7).  The  Initial  conditions  in  cylindrical 
coordinates  (r,Q)  are:  p«  1.67,  p-1,  *>0  for  r  *  l;  p*6.67,  p«4, 
m»  0  for  r  >  1.  The  gas  is  ideal  with  Y“  5/3.  The  materiel  was  divided 
into  equally  spneed  zones  of  width  .02  cm. 
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Figures  Aa  and  Ab  summarize  the  results  of  the  calculations  using 
PISCES  IDE  and  PIECES  1DL.  The  Lagrange  calculation  used  a  quadratic 
artificial  viscosity  coefficient  of  2.0  and  no  linear  artificial  vis¬ 
cosity.  The  Euler  calculation  used  a  coefficient  of  2.0  for  artificial 
viscosity.  The  problem  required  258  cycles  in  IDE  and  A92  cycles  in  1DL 
to  run  out  to  time  l.A. 

Test  Problem  5;  Elastic-Plastic  Cylindrical  Plate  Impact 

At  tine  t«0  a  steel  spherical  shell  of  thickness  37.5  (ISO  tones) 
moving  inward  with  a  radial  velocity  of  .00A  impacts  a  steel  sphere  of 
radius  12.5  (50  zones). 

Stress  profiles  from  the  IDE  and  the  1DL  calculation  are  illustrated 
in  Figures  5a  and  5b,  PISCES  1DL  used  2.0  for  a  quadratic  viscosity  coef¬ 
ficient  and  0.0  for  linear  viscosity.  PISCES  IDE  used  a  coefficient  of 
viscosity  of  2.0.  It  took  198  cycles  for  IDE  to  run  out  to  time  t »  60 
and  198  cycles  for  lDL. 

Test  Problem  6:  Spherical  Inploaion 

This  problem  is  identical  to  Test  Problem  A  except  that  the  symmetry 
is  spherical  rather  chan  cylindrical.  The  results  are  illustrated  in 
Figures  6a  and  6b. 

Test  Problem  7:  Elastic-Plastic  Spherical  Plate  Impact 

This  problem  is  identical  to  Test  Problem  5  except  that  the  symmetry 
is  spherical  rather  than  cylindrical.  The  results  are  illustrated  in 
Figures  7a  and  7b. 
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EULER  LAGRANGE  -  FLUME  TEST  B  -  LAGRANGE 
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SECTION  5 

CONCLUDING  REMARKS 

1.  The  PISCES  IDE  and  the  PISCES  1DL  code*  have  been  shown  to  be 
of  comparable  accuracy  when  each  uses  the  sane  number  of  zones.  However, 
since  for  many  problems  over  twice  os  ,mr»y  cycles  are  used  by  PISCES  1DL 
to  run  to  the  saate  problem  tine,  ut>r/;  PISCES  IDE  is  often  much  more  cost 
effective  than  using  PISCES  1DL. 

2.  The  PISCES  IDE  code  could  and  should  be  extended  in  the  near 
future  to  accurately  handle  more  than  one  material  in  a  given  problem. 

The  material  interface  must  be  treated  in  a  second-order  manner  consistent 
with  the  accuracy  of  the  rest  of  the  calculation.  Once  this  is  accosplished , 
the  PISCES  IDE  code  will  have  virtually  all  of  the  capabilities  of  the 
PISCES  1DL  code  and  may  be  used  to  solve  a  wide  class  of  problems  that 

have  traditionally  been  considered  to  be  strictly  Lagrangian  in  nature. 

3.  Once  the  smterial  interface  problem  has  been  solved  for  one 
dimension,  serious  thought  should  be  given  toward  extending  this  compu* 
tational  method  to  two  and  three  dimensions.  The  accuracy  of  the  elastic- 
plastic  calculation  combined  with  the  inherent  ability  of  an  Eulerlan 
code  to  handle  problems  Involving  large  deformations  should  make  this 
method  an  extremely  useful  computational  tool  for  multi-dimensional  calcu¬ 
lations.  It  should  be  able  to  accurately  compute  phenomena  which  have 
heretofore  been  amenable  only  to  a  coupled  Eulerlan- Lagrangian  computer 
code. 


31 


1.  Wilkins,  N.  L. ,  "Calculation  of  Elastic-Plastic  Plow,"  UCRL-7322, 

Rev.  1,  January,  1969. 

2.  MacCoruack,  R.  W. ,  "Hie  Effect  of  Viscosity  in  Hypervelocity  Iapact 
Cratering,"  AXAA  Paper  No.  69-3S4,  May,  1969. 

3.  Jardin,  S.  C. ,  "A  Numerical  Investigation  of  Unsteady  Magnetogae- 
dynsmics,"  S.H,  Thesis,  Massachusetts  Institute  of  Technology,  June 
1973. 

4.  Roache,  P.  J. ,  Computational  Fluid  Dynamics,  Hermosa  Publishers, 
Albuquerque,  N.M.  (1972). 

5.  Richtmyer,  R.  D.  end  Morton,  K.  W. ,  Difference  'Methods  for  Initial 
Value  Problems,  2nd  ed..  Interscience  Publishers ,  New  York  (1967). 

6.  Payne,  R.  B. ,  "A  Numerical  Method  for  a  Converging  Cylindrical  Sh'jck," 
J.  of  Fluid  Mechanics. (2)p.  185  (1957). 

7.  Lapldus,  A.,  "Computation  of  Radially  Symmetric  Shocked  Flows," 

J.  of  Coaqiutstional  Physics, (8) ,1,  p.  106  (1971) 


8-  ” 


.I..---. -■■v. 


*nxif**x£aujatmM;tw*jaM!xaBwmt* 


TCAM  73-14 


APPENDIX  A 

ACCURACY  AND  STABILITY  OK  THE  DIFFERENCE  EQUATIONS 


The  accuracy  and  stability  of  finite-difference  methods  such  as  those 
previously  described  cannot  be  completely  analyzed  in  the  general  non¬ 
linear  form.  The  method  of  analysis  we  will  use  is  to  first  linearize  the 
set  of  differential  equations  (1)  (with  the  vector  R  absent)  and  then  to 
study  the  amplification  of  Fourier  components  of  the  solution  by  the  dif¬ 
ference  method  applied  to  the  linearized  set. 

The  set  of  equations  we  will  analyze  is 


a  -  sa  -  =j a  - 

StU+AfcU  -  PUB^U 


where  A  and  B  are  the  Jacokisn  matrices  of  F  and  G  with  respect  to 
U  and  are  considered  to  be  constant. 
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Application  of  the  difference  method  to  the  set  (A-l)  gives 
r;  *  —  n  At^f-'n  -jnl  At;jr-*n  ;»nl 

uj  =l'i  'srA[Vx'ujJ  +  p^B|V*yjJ 
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Next,  rewrite,  letting  C  -»  A-ppB 
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and  convert  to  a  one-step  difference  method  by  substitution, 
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Truncation  error 


The  truncation  error  of  the  difference  approximation  (A-2)  to  the 
equation  (A-l)  is  readily  obtained  from  Taylor's  series  expansions.  Thus, 
letting 


»n+1  T?n  ^  ^  (At f  b2 U  , 

Uj  -uj  +  Atat+V"i?+ *•* 


n  -n  .  A  SU  .  (Axf  b2V  . 
=  u,  +  Ax-c-  +  -v wJ~  rrs  + 


3k 


and  using  the  differential  equation 


_  ,  P  au 
at  c  ax 


au 


to  eliminate  highei  order  time  derivatives  in  favor  of  spatial  derivatives, 
the  difference  equation  (A-2)  can  be  shown  to  be  equivalent  to  the  equation 


au  .  2  au  =  =* 
at +  ax  ta 


where  the  leading  terms  of  the  truncation  error  T^  are 


V-lH  £?- 


a^u  .  c 2At  r, .  ■? 

j  (Ax)  -  C  (At) 


J 


+  •« 


It  is  interesting  to  note  that  this  method  would  be  exact  (i.e.,  zero 
a  Av  =»  =» 

truncation  error)  if  C  *  -jj-I  where  I  is  the  identity  matrix.  The 
leading  tern  in  the  truncation  error  is  seen  to  be  responsible  for  dis¬ 
persion  while  the  second  term  causes  dissipation. 
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Stability 

We  investigate  the  effect  of  the  difference  equation  (A-2)  on  a  single 
Fourier  component  of  the  solution  by  making  the  substitution 


^  n  .  ~  n  ikj  Ax 
J.  ■'Dr  e  J 


(A-3) 


where  U  Is  the  amplitude  vector  of  this  Fourier  component,  i  =  <7-1  , 
k  is  the  wave  number,  and  rn*e*U,n^t  1*  the  growth  factor.  The  stability 
criterion  will  be  that  |r]  1.  If  this  condition  is  satifled,  no  har¬ 

monic  is  amplified  at  all,  whereas  if  it  is  violated,  there  is  some  har¬ 
monic  that  is  amplified  without  limit  as  n  Increases. 


Making  the  substitution  (A-3)  into  (A-2)  yields,  after  simplification, 


r  U  = 


St  A.  =*  At3  =*  2  r  1 

I  -  i  ^  Cs in(kAx)  +  C  jcos (kAx)  -  lj 


or,  upon  defining  the  "amplification  matrix"  Z  by 
—  3  -!  At  / At  32 
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Z  =  T- iC  —  Jin(kAx)  +  C  [cos(kAx)  -  l] 


we  have 


{r  I  -  Z}  •  U  -  0 

A  unique  solution  will  exist  only  if  the  determinant  of  the  matrix  in 
brackets  vanishes,  i.e., 


det{r I  -  z} 
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This  is  equivalent  to  the  stetenent  that  r  must  equal  one  of  the  eigen- 
s*  rt 

values  of  Che  matrix  Z.  The  eigenvalues  of  the  matrix  Z  may  be  expres¬ 
sed  in  terms  of  the  eigenvalues  of  the  matrix  C  by  making  use  of  the 
Spectral  Mapping  Theorem  which  states  that  if  Z  »  /(C)  is  a  rational 
function  of  the  matrix  C  and  are  the  eigenvalues  of  C,  then  /(Xc) 
are  the  eigenvalues  of  Z.  Thus, 

*  «  1  '  i  sin(kAx)  +  X*[coa(kAx)  -  l] 

taking  the  square  of  the  absolute  magnitude  and  Imposing  the  condition 
|  r  |  £  1  gives  the  requirement 


2  /  X  At,2  /  A  OCV- 

|r|  «  1+(“3r)  Cco.<k/ix)  -  13  +  1-^)  sin2(kAx)£l 


This  is  seen  to  trace  out  an  ellipse  in  the  complex  plane  as  (kAx)  varies 
from  0  to  2n. 
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It  is  easiest  to  evaluate  this  determinant  by  going  back  to  the 
primitive  variables  p,  u-m /p,  e »  (e/p)  -  u2/2,  and  s^^s/p.  Making 
use  of  the  relations 


&  .  &  de  „  1  & 
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Equation  (A-6)  can  be  factored  to  give 


(Xc-u>3[Xc-2uXc  +  u3-Ip  -  !?"  -  0 

=t 

Thus  three  of  the  eigenvalues  of  the  matrix  C  are  degenerate  and  are 
equal  to  the  particle  velocity  u.  The  other  two  are  obtained  by  solving 
the  above  quadratic.  This  gives 


u 

u 

u 

u  +  c 
u  -  c 


(*-7) 


where 


(A-8) 


The  quantity  c  is  seen  to  be  the  sonic  velocity  for  a  general  equation 
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APPENDIX  B 
INPUT 


The  PISCES  IDE  code  is  input  exactly  the  cine  es  the  PISCES  1DL  code 
with  the  following  exceptions  and/or  restrictions. 

1.  A  single  TYPE  57  card  aust  be  input  with  the  number  1.0  in 
the  first  field.  This  tells  the  code  to  do  an  Eulerien 
rather  then  a  Lagrangian  calculation. 

2.  Only  one  awterial  type  is  allowed  in  a  problen  and  it  aust 
have  material  number  1. 

3.  The  Boundary  Options  card  (TYPE  41)  must  be  Included  in  the 
input  but  it  is  ignored  by  the  code.  Standard  boundary  con¬ 
ditions,  unless  otherwise  supplied  are: 

LEFT  BOUNDARY:  Kept  at  initial  conditions  for  plane 
and  flume  symmetries,  symmetry  axis 
for  cylindrical  and  spherical  sym¬ 
metries. 

RIGHT  BOUNDARY:  Zero  gradients  at  grid  point  (JMAX  -  1) 

4.  Yield  models  based  on  the  distortional  energy  may  not  be  used 
as  the  distortional  energy  is  not  explicitly  calculated  by  the 
code. 

5.  The  shear  modulus  is  assumed  to  be  constant  and  to  be  equal 
to  GO. 

6.  The  input  parameter  DISL  must  always  be  included,  but  it  is 
Ignored  for  cylindrical  or  spherical  symmetry  problems. 
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Instead,  DISL  Is  computed  by  the  code  so  that  the  first  two 
grid  points  are  centered  about  the  symmetry  point.  The 
values  at  the  first  grid  point  are  then  not  to  be  interpreted 
as  being  physical,  but  are  chosen  to  provide  a  "symmetry" 
boundary  condition. 

7.  The  code  has  perhaps  its  greatest  usefulness  in  calculating 
"moving  wall  flume"  (NSWD*5)  problems.  In  these  calcula¬ 
tions,  Subroutine  MOVWAL  is  not  called  as  it  is  in  PISCES 
1DL.  Instead,  the  user  must  supply  a  Subroutine  AREE(AT,AX, 
AI,DH)  which  is  called  twice  for  each  zone  for  each  cycle. 

The  variables  AI  and  DH  are  not  presently  used  by  the  code. 
(These  should  be  set  to  1.0  in  the  subroutine.)  The  variable 
AT  must  return  the  derivative  of  the  area  at  zone  J  with  res¬ 
pect  to  the  time,  divided  by  the  area.  Similarly,  the  vari¬ 
able  AX  must  return  the  derivative  of  the  area  at  zone  J  with 
respect  to  r,  divided  by  the  area.  For  the  method  to  by  ful¬ 
ly  second  order,  these  derivatives  must  be  evaluated  at  time 
n  on  the  first  call  and  time  n+  1  on  the  second  call. 


5-  \o 


I 


5 


t 


\ 


APPENDIX  B 

PISCES™  1DL 
MANUAL  A 


l 


.K 

i 


b-  H 


.Cx-W-'  V.  C*v  - 


^»««s«»sov«asiB3S5Si®g  r 


1 1)1. -A 


CONTENTS 


GENERAL  DESCRIPTION  I 

Equations  of  Motion  I 

Constitutive  Equations  3 

Boundary  and  Initial  Conditions  4 

Void  Option  5 

S|>ccial  File  Option  5 

Flume  Option  5 

Static  Option  5 

Extending  Re/one  Option  5 

Microzone  Option  6 

H«.-at  Flow  Option  6 

Restart  Option  A 

Data  Display  Option  6 

FINITE-DIFFERENCE  EQUATIONS  7 

APPENDIX  A:  Derivation,  of  Equation  of  Motion  in 

Flume  Symmetry  23 

APPENDIX  B:  Relations  Between  Measures  of  Volume  Change  25 

APPENDIX  C:  The  q  Stability  Criterion  27 

APPENDIX  D:  Accuracy  of  the  Pressure-Energy  Iteration  31 

APPENDIX  E:  Relation  Between  er  ,  e,  and  Compression 

in  Plane,  Cylindrical  and  Spherical  Geometry  33 

APPENDIX  F:  Formula*,  for  Sound  Speed  in  Elastic-Plastic  Material  37 


iii 

6. I* 


iSfjJfcrTl'r 


I  DIM 


GENERAL  DESCRIPTION 


EQUATIONS  01'  MOTION 


PISCES  IDL  is  a  one-dimensional.  time-dependent  finite-difference  Lagrange  aide 
which  is  used  to  calculate  the  dynamic  motion  of  contmoa  including  elastic-plastic  and 
hydrodynamic  media.  The  differentia!  equations  that  govern  the  motions  of  these  media 
are  approximated  by  difference  equations  applied  to  a  network  of  /ones  that  describe 
the  physical  space  occupied  by  the  media.  Groups  of  contiguous  /.ones  of  the  same  material 
are  described  by  constitutive  equations  which  are  coupled  to  the  equations  of  motion. 
Time  is  also  discretized  and  its  formulation  insures  stable  difference  equations.  The 
Lagrangian  formulation  of  the  differentia!  equations  requires  that  each  zone  contain  a 
constant  mass  element  of  material  that  moves  and  distorts  in  space  and  time,  the  motion 
for  ali  the  zones  approximating  the  continuum  motion.  In  a  one-dimensional  formulation 
only  one  space  variable  exists,  symmetry  considerations  account  for  the  other  two.  There 
are  three  possible  symmetries  plane,  cylindrical,  and  spherical,  in  plane  symmetry,  zones 
are  infinite  slabs  and  the  thickness  of  a  slab  is  the  ?one  si/e.  In  cylindrical  symmetry, 
/ones  are  concentric  hollow  cylinders  of  infinite  length.  Spherically  symmetric  zones  are 
concentric  hollow  spheres.  In  both  of  these  latter  geometries,  zonal  widths  are  radial 
thicknesses  of  the  cylinders  and  spheres.  res|iectively.  The  partial  differential  equations 
that  govern  the  motions  in  PISCES  IDL  are  shown  on  the  following  pages. 
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Continuity  Equations 
Plane  symmetry 
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Energy  Equations 


Plane  symmetry 


E  -  -<P  ♦  q)V  +  v/*r  4  H 


Cylindrical  symmetry 


Spherical  symmetry 

E  -  -(p  +  q)V  +  Ji  +  2sc  +  II 
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r  stress  devMim 


t  stress  dcviat'ii 
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p  =  hydrostatic  pressure 

q  =  artificial  viscosity 

_  po 

V  =  relative  volume  =  —  ,  pQ  =  reference  density 

E  =  internal  energy  per  original  volume,  H  =  heat  flow. 


In  the  equations  of  motion  the  position  and  velocity  of  a  particular  zone  arc  defined 
at  the  zone  corners  (grid  points  that  define  the  /.on;;)  boundary)  while  other  variables 
such  as  pressure,  stress,  energy,  etc.,  are  defined  as  averages  over  the  zonal  interior.  These 
differential  equations  are  approximated  by  the  finite-difference  equations  developed  by 
Wilkins  (Reference  I). 

CONSTITUTIVE  EQUATIONS 

Besides  the  equations  of  motion  it  is  necessary  to  have  constitutive  equations  that 
describe  the  media.  These  material  description  equations  are  the  equations  of  state, 
stress-strain  models  and  transport  property  descriptions.  In  general,  these  functional 
relationships  may  take  any  form-discontinuous  functions  as  well  as  tables  are  used  to 
describe  phase  changes,  hysteresis,  cracking,  etc. 

The  equations  of  state  calculate  pressure  from  the  density  and  the  total  internal 
energy  density  of  a  zone  in  a  two-step  iteration.  The  interim  internal  energy  density  is 
evaluated  from  the  old  internal  energy  density,  the  new  heat  flux,  the  old  pressure,  and 
the  new  change  of  volume  in  the  zone.  This  first  approximation  of  the  internal  energy 
density  is  then  used  in  the  equation  of  state  to  get  a  first  approximation  of  the  new 
pressure.  A  revised  estimate  of  the  internal  energy  density  comes  from  tms  first 
approximation  of  the  hew  pressure  and  the  new  distortional  energy  density.  Finally,  the 
new  pressure  in  a' zone  is  determined  from  the  second  approximation  to  the  internal  energy 
density.  In  this  way  a  simultaneous  solution  of  the  constitutive  relations  and  the  second 
law  of  thermodynamics  is  effected  independently  of  the  form  of  the  equation  of  state. 

An  example  of  an  available  stress-strain  model  is  the  so-called  elastic-plastic  yield 
stress  model.  A  yield  stress  model  is  a  law  that  regulates  the  stress  deviators.  In  a  calculation, 
the  stresses  are  decomposed  into  a  hydrostatic  component  and  a  deviatoric  component. 
The  stress  deviator  describes  the  resistance  of  the  material  tc  shear  distortion  and  is 
calculated  in  terms  of  an  incremental  stress  resulting  from  an  incremental  strain.  The 
deviators  are  limited  by  the  yield  stress  according  to  some  rule  that  describes  the 
elastic-plastic  behavior  of  the  material,  in  this  model,  the  deviatoric  behavior  of  materials 
is  determined  by  the  shear  modulus  and  the  yield  stress.  (For  a  fluid,  both  are  zero.) 
A  icmperature-dependent  sheer  modulus  is  available  and  should  be  chosen  to  be  consistent 
with  the  bulk  modulus  and  Poisson's  ratio  for  the  material.  The  yield  stress  is  defined 
as  the  stress  at  Ihc  yielding  point  in  a  uniaxial  stress  test  and  may  he  dependent  on 
thermodynamic  variables.  For  example,  the  code  has  available  work-hardening  models  in 
which  yield  stress  is  a  function  of  distortional  energy.  Mohr-Couloinh  models  in  which 
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Time-dependent  pressure  profiles  as  well  as  time-dependent  velocity  profiles  (pistons) 
can  be  automatically  applied  to  either  boundary.  Fixed  and  free  boundary  conditions  are 
merely  sub-cases  of  the  above  capabilities.  Energy  boundary  conditions  are  specified  by 
applying  temperature  and  heat  flux  histories  at  the  appropriate  boundaries.  Internal  or 
kinetic  energy  can  be  deposited  into  any  region  as  an  initial  condition.  Time-  and 
space-dependent  energy  deposition  and  a  zone-by-zone  detonation  description  may  be 
included  in  the  solution  of  a  problem  as  required.  The  code  logic  simulates  a 
Chapman-Jouguet  detonation. 

Once  the  geometry,  material  properties,  initial  conditions,  and  houndary 
conditions  are  specified,  the  calculation  of  the  problem  is  ready  to  proceed 
Tin  first  lime  step  that  the  problem  takes  must  be  specified,  but  Irons  that 
point  on.  the  code  automatically  calculates  the  subsequent  time  steps  accord¬ 
ing  to  stability  cntoiia  imposed  bv  tnc  linite-differenee  equations.  Physically, 
the  restriction  imposed  on  the  calculation  is  that  no  signal  may  cross  a  zone 
in  one  time  step  so  that  a  zone  reacts  only  to  activity  being  transmitted  from 
a  neighboring  zone. 

The  accuracy  of  the  solution  in  a  i.agrangian  code  depends  on  the 
density  of  zones.  With  coarse  zoning,  gross  features  of  the  solution  such  as 
total  energy  and  momentum  in  a  region  may  be  realistic,  but  the  actual  values 
of  a  variable  at  a  point  are  apt  to  be  in  error.  The  resolution  in  both  space  and 
time  is  proportional  to  the  zone  size.  A  strong  shock  is  handled  by  the  finitc- 
ditference  equations  by  using  Von  Neumann's  artificial  quadratic  viscosity 
(Reference  2).  which  is  introduced  to  spread  any  shock  front  over  a  distance 
of  about  three  zones.  Zone-to-zone  noise  may  be  smoothed  by  an  artificial 
linear  viscosity  logic. 
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PISCES  IDL  contains  several  options  that  are  worth  special  mention.  Among  these  | 

are  the  void,  remap,  flume,  static,  extending  rezonc,  microzone,  heat  llow.  restart,  and  1 

data  display  options.  These  capabilities  are  described  as  follows:  | 

VOID  OPTION  I 

Between  any  two  mass  regions  it  is  possible  to  define  a  void  of  finite  thickness. 

Motion  in  these  regions  can  expand  or  shrink  the  void.  However,  once  Hi'.’  void  is  closed, 
it  cannot  he  reopened.  This  option  is  valuable  when  dealing  with  particle  interaction 
problems  and  would  probably  be  useful  when  dealing  with  some  types  of  structural 
response. 

SPECIAL  FILE  OPTION 

The  special  file  capability  allows  information  calculated  by  PISCES  IDL  to  be 
punched  onto  computer  cards,  written  onto  tape  or  disk,  or  stored  by  any  computer  re¬ 
source  on  its  media.  This  output  can  be  formated  in  such  a  way  as  to  be  input  to  another 
code. 

FLUME  OPTION 

The  flume  option  turns  the  code  into  a  one-and-one-half-dimensional  code  by  defining 
a  variable  area  duet  as  a  boundary  condition  along  the  direction  of  How.  The  duet  may  t 

be  fixed  or  free  to  respond  to  the  flow  pressures.  In  the  latter  case  mot  ion  of  the  duct  < 

walls  may  be  programmed  to  have  some  pre-described  resistance  or  motion.  A  real  strength  j 

wall  model  is  available  upon  request.  For  material  moving  in  a  direction  parallel  to  a 
fixed  or  slowly  moving  boundary,  the  Hume  option  may  obviate  the  need  for  a  more  I 

expensive  two-dimensional  calculation.  Flume  materials  must  always  he  Hinds  while  wall 
materials  may  be  solid.  j 

STATIC  OPTION 


Nondynamic,  equilibrium  stress  distributions  may  be  computed  using  PISCES  IDL. 
An  input  switch  changes  the  time  variable  to  an  iteration  variable  converting  the  code 
to  a  stre.>s  diffusion  equation  solver  rather  than  a  stress  wave  equation  code.  This  option 
is  relatively  new  but  very  powerful. 

EXTENDING  REZONF.  OPTION 

The  rezoner  allows  a  user  to  extend  automatically  the  spatial  domain  included  in 
a  calculation  while  a  problem  is  in  progress.  For  explosions  emanating  from  a  point,  line, 
or  plane,  the  spatial  length  of  the  problem  is  doubled  when  activity  nears  the  outside 
boundary.  The  total  number  of  zones  with  which  the  problem  starts  is  liekl  fixed  so 
that  the  effect  of  this  rezone  procedure  is  that  all  zones  are  approximately  doubled  in 
size.  There  is  no  limit  to  the  number  of  rezonings  that  may  be  allowed  in  a  problem. 
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bul  since  l he  accuracy  of  a  solution  depends  on  the  zone  size,  resolution  can  become 
poor  alter  several  rc/.onings.  Generally,  for  one-dimensional  explosions,  this  procedure  lias 
the  advantage  of  having  the  finest  zoning  when  it  is  most  necessary. 

MICROZONE  OPTION 

The  micro/.oner  is  used  when  fine  zoning  requirements  at  a  shock  front  are 
economically  too  stringent  for  the  entire  grid.  This  option  allows  fine  zoning  to  be  initially 
specified  only  where  it  is  necessary,  while  coarse  zoning  is  defined  elsewhere.  Then  the 
code  automatically  moves  the  microzoned  (finely  zoned)  region  to  follow  the  activity 
that  needs  the  fine  zoning.  A  single  pulse  in  a  slab  where  the  entire  pulse  is  small  compared 
to  the  width  of  the  slab  is  an  example  of  a  problem  that  may  need  the  micro/oning 
option. 

HEAT  FLOW  OPTION 

Diffusion  logic  is  available  to  handle  heat  conduction  problems  l<  mi  id's  law  is 
coupled  to  the  difference  equations  of  material  motion  under  the  assumption  that  the 
heat  conduction  need  not  be  time  centered  in  the  energy  equation  when  momentum  el  feels 
control  the  time  step. 

RESTART  OPTION 

The  restart  option  allows  the  user  to  complete  a  problem  in  sequential  steps.  Edits, 
rezones,  boundary  conditions,  and  material  parameters  may  all  be  redefined  at  restart  time. 

DATA  DISPLAY  OPTION 

Automatic  piot  display  is  available  in  PISCES  IDL  through  the  PI-PLOT1  M  Plotting 
System.  The  plotting  capability  allows  parametric  plots  of  any  variable  as  a  function  of 
any  variable  at  a  point  in  Lagrange  space  with  time  as  a  parameter  or  at  a  particular 
time  with  laboratory  coordinate  as  a  parameter.  These  plots  may  be  plotted  automatically 
on  a  Cal  Comp  incremental  or  zip-mode  plotter,  an  SC-4020  CRT  plotter,  or  on  the  printer, 
in  addition,  shock  front  plots  (values  of  a  variable  at  the  shock  front  as  a  function  of 
distance  and  time)  are  available  and  are  plotted  on  a  log-log  scale  for  convenience.  Peak 
value  printouts  and  plots  are  also  available. 

Automatic  print  edits  arc  available  as  a  function  of  problem  time,  cycle  number 
and  activity  front 
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FINITE-DIFFERENCE  EQUATIONS 


Thu  differentia!  equations  solved  in  PISCES  IDL  have  been  summarized  in  the  general 
description.  The  first  order  explicit  finite-difference  equations  that  have  been  used  to 
approximate  the  differential  equations  are  discussed  below. 

For  the  following  discussion  there  are  a  few  definitions  and  conventions  which  must 
be  addressed.  First.  PISCES  IDL  is  a  one-dimensional  Lagrange  code.  The  space  variable 
is  divided  into  /.ones  of  fixed  mass.  Grid  points  bound  each  /.one.  An  interior  grid  point 
is  bounded  by  a  /.one  on  either  side  while  a  boundary  grid  point  has  one  of  the  /ones 
replaced  by  a  boundary  condition.  A  left  boundary  grid  point  has  a  /.one  only  on  its 
right:  a  right  boundary  grid  point  has  a  /.one  only  on  its  left.  Indexing  of  grid  points 
is  generally  denoted  by  j  and  increases  from  left  to  right.  Indexing  of  zones  is  done  by 
associating  a  zone  with  its  right  side  grid  point  index.  Thus,  the  leftmost  boundary  grid 
point  (or  the  lowest  index  number)  has  no  /.one  associated  with  it. 

The  calculations  are  divided  into  two  types-grid  point  and  zone.  Associated  with 
each  interior  grid  point  is  a  Lagrange  position,  velocity,  acceleration,  half  of  the  mass 
of  each  of  the  neighboring  zones,  and  a  heat  flow  area.  Associated  with  each  zone  is 
the  fixed  zone  mass,  relative  volume,  density,  compression,  pressure,  artificial  viscosity, 
stress  deviators,  yield  stress,  internal  energy  density,  distortional  energy  density,  sound 
speed,  and  temperature. 

Two  basic  conservation  equations  are  solved  one  for  each  of  the  two  types  of  calcula¬ 
tions.  The  Lagrange  equation  of  motion,  derived  from  conservation  of  momentum  and 
conservation  of  mass,  is  used  to  update  grid  point  variables.  The  Ugrange  thermal  energy 
equation,  derived  from  conservation  of  energy  and  continuity,  is  used  to  update  /one 
variables.* 

Updating  of  grid  point  variables  is  done  in  subroutine  MOTION.  First  the  equation 
of  motion  is  used  to  calculate  grid  point  accelerations  of  the  previous  cycle  from  stresses 
calculated  in  the  previous  cycle. 


For  NSWD  *  I  (plane  symmetry) 
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For  NSWD  *  2  (cylindrical  symmetry) 
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*  Deri  vat  inn  of  basic  equations  is  available  upon  request. 
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For  NSWD  =  3  (spherical  symmetry) 
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For  NSWD  *  4.5  (flume  svmmetrv)* 
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A  general  form  of  these  c«|uations  may  he  written  in  the  following  way: 
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where  d  *  NSWD  and  for  NSWD  >  3.  set  d  -  I.  This  equation  may  be  rewritten  in 
terms  of  the  relative  volume.  vn  s  p  /p”  where  pQ  is  a  reference  density  usually  set 
equal  to  the  initial  density  0 
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The  finite  difference  analog  of  the  coefficient  of  dZ^  is  the  reciprocal  of  the  grid  point 
areal  density.  The  areal  density  to  the  lei t  of  the  grid  point  j  and  the  areal  density 
to  the  right  of  the  grid  point  j  are.  respectively. 
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*Sve  Appendix  A  Jar  derivation. 
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The  stress  difference  in  brackets  is  calculated  from  an  average  of  the  stress  difference  in  the 
zones  to  the  left  and  right  of  the  grid  point.  These  are.  respectively. 
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Thus,  for  an  interior  grid  point,  i.e..  a  non-boundary  grid  point,  the  finite  difference  formula 
for  the  equation  of  motion  is 
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When  calculating  a  free  left  boundary,  m”  is  zero  and  E91  is. zero.  For  a  pressure 
time  history  applied  to  the  left  boundary,  is  still  zero  but  takes  on  the  values 
of  the  prescribed  pressure.  For  a  fixed  boundary,  the  above  equation  is  superfluous  and 
r j+^  is  set  to  zero.  Similarly,  a  velocity  condition  on  the  boundary  (piston)  is 
automatically  set  equal  to  r?*?  Right  boundaries  are  handled  in  a  similar  way,  but 
m”  and  X^2  are  adjusted.  'Only  one  difference  exists,  J^2  takes  on  values  of  minus 
the  pressure  time  history  when  a  pressure  profile  is  used. 

Void  boundaries  are  calculated  us  free  boundaries.  A  left  void  boundary  is  a  right 
free  boundary  of  the  solid  region  on  the  left.  A  right  void  boundary  is  a  left  free  boundary 
of  the  solid  region  on  the  right.  Void  closure  results  in  the  two  boundary  points  becoming 
one  interior  point.  The  velocity  of  closure  is  computed  to  give  exact  closure  while 
conserving  momentum.  This  will  result  in  an  induced  energy  error.  To  see  this,  consider 
the  formulas  for  closure  velocity  derived  from  conservation  of  momentum 
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and  from  conservation  of  energy 
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the  formulas  become 
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Furthermore,  if  r^  -  0.  then 
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The  differences  are  clear. 


When  Hi|uation  (I)  is  used  to  get  acceleration.  r"j.  subsequent  calculations  for  new 
velocity.  f"+^.  and  new  position.  r"+l.  are  computed  from  the  following  time-centered 


equations: 

•  r.+lj  .n-ls  .  -n.^n 
r  ’  =*  r  ^  +  r  At 


n+i  n  ,  •n+lj.#.n+l( 
rj  “  rj +  rj  “  • 


Once  the  positions  and  velocities  of  grid  points  forming  a  /one  have  been  updated 
in  subroutine  MOTION  (and/or  VOID),  pressure,  volume,  energy,  etc.,  can  be  calculated 
for  the  zone.  These  and  other  zone  interior  variables  are  computed  in  subroutine  /.ONI-.. 


The  first  zone  interior  variable  calculated  in  subroutine  /ONI:  is  the  relative  volume. 
It  should  be  obvious  that  once  new  coordinates  are  available,  one  can  compute  a  new 
volume. 


For  NSWD  *  I  (plane  symmetry) 


TRUE  VOLUME"*1  -  r^-r"41 
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In  cases  where  a  velocity  boundary  condition  is  used,  only  bquation  (3)  applies.  Tltcnr"  ^ 
and  rj  are  the  velocity  and  position  from  the  previous  cycle  (or  initial  conditions);  At"  and 
Atn*^  come  from  stability  considerations  based  on  the  previous  cycle  (or  an  initial  value). 


For  NSWD  *  2  (cylindrical  symmetry) 


TRUE  VOLUME 
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For  NSWD  *  3  (spherical  symmetry) 


TRUE  VOLUME^1  -  |n  j  -  (r£$]< 


For  NSWD  *  4  (fixed  wall  Hume)* 


J 

TRUE  VOLUME?* 1  -  £  (frustums  defined  by  Rj.) 


where  R^  are  wall  points  between  r*jtj  andrj*1, 
For  NSWD  *  5  (Moving  wall  flume)* 


J 

TRUE  VOLUME"*1  -  J  (frustums  defined  by  R^1) 


where  R?*1  are  wall  points  between  r?*1  and  r?* 1 . 
i  j-i  j 

To  compute  relative  volume,  V .  form  the  expression 


[TRUE  VOLUME j  j[ REFERENCE  DENSITY J 
i  MASS ^ 

Vice  Appendix  A  for  a  picture  of  a  typical  flume  .section. 
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Computing  new  relative  volume.  V0**.  ‘•’O'*;  t*’.e  change  i'l'  relative  volume.  AVnVi. 
itnild  also  introduce  round-off  error  if  ’’  <'  V1*41.  In  general.  Vn+1  is  near  1.0 

and  h!n*'i  is  very  small.  To  elunira'.  Me  possibility  ol  this  kind  of  rotind-olf.  a 
transformation  of  variables  is  made,  "’’he  variable  compression.  C.  is  defined  as  follows:* 


C  -  1.  -V 


AC  =  A(l.  -  V)  .  -AV.  (4) 

Compression  and  change  in  compression  are  generally  both  numbers  of  the  same  order 
of  magnitude.  The  equation 


therefore,  does  not  introduce  a  round-off  error.  New  relative  volume  comes  directly  from 

v”*1  «  1  -c**1 

vj  1.  Cj  . 

The  formula  lor  change  in  compression  in  plane  symmetry  is 


What  we  have  here  is  a  formula  for  the  change  of  relative  volume  that  depends  on  velocities, 
not  position.  If  the  velocities  yield  no  change  in  relative  volume,  then  one  is  guaranteed 
that  the  distance  between  zones  actually  is  unchanged  since  new  positions  are  computed 
from  new  velocities. 

To  use  this  formula,  one  must  be  sure  that  it  agrees  with  the  development  based 
on  coordinates.  This  must  be  a  requirement  since  the  TRUE  VOLUME  formulas  are  exact 
for  large  displacements.  Again,  looking  at  the  plane  symmetry  case,  we  sec 

p.  /_mi  n\  /  n+’  n  v 

m 

which  reduces  to  the  result  obtained  previously. 

*See  Appendix  B  for  relations  between  volume,  relative  volume,  dilation  and  compression. 


£.13 


f 


Finally,  in  spherical  geometry  the  change  of  relative  volume  from  coordinates  is 

“r4-^3"|[('rv(vf)-[H)'-(^)]| 

<-^^|l(^-ra1-|('r-i)°-«-.)]|- 


From  velocities 


K) 

j 

d 

lin 

|/rn+1\3  . 

/rn  fl 

dt 

■i 

dt 

1  3 

[\rj  ) 

\  J-V  ] 

-r*-^*-|,T(*r4f  '«ss«3/!  *■■ 


If  this  formula  is  expanded,  it  will  not  agree  with  the  formula  derived  directly  from 
coordinates.  There  is  an  error  reflecting  the  lack  of  second-order  accuracy  in  taking  the 
derivative,  lhe  formula  in  terms  of  velocities  becomes 


•  #  (tf)’  *  |  j,°*4 


^  \('7hj  -  (^)1 


when  the  error  is  accounted  for.  The  change  in  compression  is 

<  ■  ^  >(»T(*r*y  •  *£(■?!)  *  *7 


d-  '5 


'if. 


I1)L-A 


where  again 


*f  -  {(<:i  -  (^f\- 

New  relative  volume  is  now  computed  from  new  compre.wion. 

cj*‘  »  Cj  +  t£^ 


..n*1  ,  „n+! 

Vj  ml-°*  * 


(5) 

(6) 


Change  in  relative  volume  is  of  course  the  negative  of  the  change  in  compression, 

AvJ+*  .  -  AcJ*%  .  (7) 

The  rate  of  change  of  relative  volume  can  be  computed  directly. 


v"** 

^5- 


(*) 


j 


j 


Once  the  new  relative  volume  has  been  determined,  the  artificial  viscosity,  pressure,  energy, 
stress  deviators  and  new  time  step  may  be  determined. 


Artificial  viscosity,  q”*  *.  is  computed  in  tin:  following  way  (Appendix  C): 


n+%  n+%  n4% 

•»  j  "  v  +  V 


(9> 


j 


The  artificial  viscosity  is  the  sum  of  an  artificial  quadratic  viscosity,  q”*^ .  and  an  artificial 
linear  viscosity,  j 

Lj 


.  pc,  ,,  _x2  v"*  ,vr* 

sSf  *  -  \(ir  )  |-«5 


II. 


(9.1) 


where  pQ  =  reference  density,  Vj+^  =  relative  volume,  L r  is  distance  sciov*  the  /.one, 
and  Cq  =  coefficient  of  quadratic  viscosity. 

P0  « 

J  vj  J  vi  1 

where  =  coefficient  of  linear  viscosity  and  c  =  sound  speed. 

The  computation  of  pressure  and  energy  cannot  be  done  in  a  straightforward  manner 
because  one  depends  on  the  other.  Therefore,  an  iterative  approach*  is  required.  First, 
a  preliminary  calculation  is  made,  estimating  the  new  energy  from  new  change  in  relative 
volume,  new  artificial  viscosity,  and  oid  pressure.  Distortional  energy  change  is  ignored 
since  new  stress  deviators  have  not  yet  been  calculated.  This  first  estimate  ol  the  new 
energy  is  called  En+1. 


-n+ 1  _n  1 1  /  n  .  n\  .  n+*j  ...nH 

•j  ■  Ej  -  [n’j +  pj) +  "j  j  4Vj 


If  energy  is  deposited  from  outside  the  grid,  it  is  added  into  En  prior  to  the  first 
approximation  to  the  energy,  e.g.. 


Ej  -  E"  +  I  E'At1** 


where  Ej  is  the  total  energy  to  be  deposited  in  a  i.one. 


Now  using  the  energy  estimate,  E”  ,  subroutine  STATF.  is  called  to  get  an  estimate 
of  pressure,  'pj* 1 , 


*7‘  -  /(if*,  v"1).. 


Returning  to  subroutine  ZONL,  stress  deviators  are  computed  by  first  calculating 
the  ‘‘clastic"  deviators  at  plus  time  from  strain  rates  at  half  time.  The  elastic 
one-dimensional  deviators  are  computed  from 


-  .2  +  2ofc+%  ^ \ 


3 

J 


* Appendix  D. 
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KT  -K'l'-M- 

This  quantity  is  compared  with  the  computed  allowable  yield  stress  squared.  (Yj+1)* 

/  n+i  _n+i  „mi  _n  mx\ 

Yj  UPj  *  Ej  *  Vj  '  Zj*  *  /  * 

(Notice  that  p"+1  and  e”+1  are  really  p?+1  and  E^+1.) 

The  comparison  is  made  in  the  following  way: 


It'  AYn+1  is  negative,  the  deviators  are  correct  already  everything  is  elastic.  If  AYj 
is  positive,  the  deviators  are  adjusted  by  the  ratio  of  allowable  to  elastic  yield  stress 


■  <*'/<' 

j  1  °j 


ml  ,,mi  ml 

8r  -  Y*  ®r 

j  J  j 


mi  „mi  mi 
st  *  Y*  8t  ' 

j  i  i 


(16) 


(If 


f  Y?+l  ■  0. ,  then  s”+1  »  a”*1  *  s”+1  =*  0.  exactly.  This  is  the  condition  for 
J  j  j  i 


hydrodynamic  flow.) 

From  the  new  stress  deviators  one  can  calculate  the  change  in  distortional  energy. 

.  V?*  Atn+*  Lm*  bn">  +  (d  -  1),“**  07) 

J  J  lrJ  rj  lJ  Cj  i 

Where 


C.  i<« 
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Newton's  second  law  lakes  Hie  form 


pn[R(r*)]  dr  r  *=  -S(r)n[R<r)  ]%  £(r  +  dr)n[R(r  +  dr) ]? 

-£(r*)n[R(r)  +  R(r +  dr) ][R(r +  dr)  -  R(r) ] 


Letting 


R(r+dr)  *  R(r)+dr-^  ,  etc.  gives 

.0R3drr  =  dr  R2  -  2dr  2r  +  2dr  Er 

or  or  Sr 


Nii|;iliun: 
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RELATIONS  BETWEEN  MEASURES  OF  VOLUME  CHANGE 


inittiil  volume 


actual  volume 


relative  volume  =  — 
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cubical  dilation  =  — —  1  -  — - — 
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APPENDIX  C 

THE  q  STABILITY  CRITERION 

The  derivation  of  the  q  stability  criterion  may  be  found  in  the  references  listed 
at  the  end  of  this  appendix.  When  reading  these  references,  it  will  help  to  know  the 
corresponding  notation.  PISCES  1DL  Cq  is  c  in  Von  Neumann  and  Richtmyer’s  article 
and  a  in  Richtmyer’s  book.  The  8  in  Richtmyer’s  book  is  Cq  times  the  zone  size  ahead 
of  the  shock.  Richtmyer’s  V  is  specific  volume,  and  his  E  is  internal  energy  per  unit 
mass.  Richtmyer’s  R  is  the  same  as  PISCES  IDLr.  but  for  the  distance  across  a  zone 
he  often  uses  the  more  complicated  expression, 


where  r  is  the  original  coordinate  and  a  is  1 ,  2,  or  3  for  plane,  cylindrical,  or  spherical 
symmetry,  respectively. 

A  simplified  derivation  of  the  quadratic  qq  stability  criterion  will  be  given  here. 
The  same  method  may  be  used  to  derive  the  linear  q^  criterion.  The  differential  equations 
in  plane  symmetry  are 


i  csv  ay 
v  at  “  ax 


where  u  is  velocity  and  Ax  is  the  zone  size.  We  want  to  examine  the  case  where  pressure 
is  negligible  compared  to  qq.  Set  p  equal  to  zero  in  the  first  equation.  Eliminate  relative 
volume  from  the  last  two  equations. 


These  equations  may  be  combined  into  one: 


,au.  ay] 
■ax1  a*  • 


Preceding  pace  Milk 
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In  order  to  make  this  equation  look  like  the  diffusion  equation,  we  want  to  take  the 
density  out  from  inside  the  derivative  on  the  right.  Say  that  the  variation  of  density  is 
small  compared  to  the  other  derivatives  (Von  Neumann  and  Riehtmyei  make  the  same 
gross  assumption), 


p4S  ■  <%U*f  r>: 


a_  au 
ax  1  ax1  ax 


^fp2|||0. 


Rearranging. 


du  „„»/*  v2  l^t  8-U 

af=2yfix)  I- 1  ^ 


Again  following  Von  Neumann  and  Richtmyer.  say  that  the  first  derivative  of  velocity 
varies  slowly  compared  to  the  second  derivative.  Then  the  last  equation  is  the  diffusion 
equation 

dU  „  d3U 

at  a*8 

with  the  diffusion  coefficient 
0  *  2C*<Ak>’  if  j 

*  2 C*«*f  !~|. 

In  Chapters  i  and  8  of  Richtmyer’s  book,  you  will  find  that  all  explicit  finite  difference 
approximations  for  the  diffusion  equation  must  satisfy  the  stability  criterion 


Applied  to  the  present  case,  the  diffusion  stability  criterion  is 


*<$!?! 


this  result  may  be  found  in  Von  Neumann  and  Richtmyer’s  article  at  liquation  <»S  with 
the  diffusion  coefficient  defined  at  Equation  45.  In  the  first  edition  of  Ricblnt>er\  l>ix>k. 
this  result  is  liquation  10.32  and  in  the  second  edition  it  is  Equation  12.48. 


:k 
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Von  Neumann  and  Richtmyer  did  not  stop  their  derivation  with  this  result,  which 
applies  in  general,  but  went  on  to  relate  this  criterion  to  the  .sound  speed  behind  the 
shock  front  for  the  case  of  a  steady  uniform  shock  in  a  gamma-law  gas.  Tlicir  final  result 
depends  on  the  shock  strength,  on  the  equation  of  state,  and  on  the  assumption  that 
the  state  is  uniform  behind  the  shock.  The  lesson  to  be  learned  here  is  that  even  the 
greatest  mathematicians  sometimes  do  not  know  when  to  stop  an  analytic  derivation  of 
a  relation  to  be  used  in  a  general  purpose  finite  difference  code.  The  simplest,  most  general 
relation  should  be  used.  If  one  works  toward  a  more  specific  analytic  result  to  put  into 
a  code,  the  code  will  be  more  complicated,  less  flexible,  and  less  general. 
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APPENDIX  D 

ACCURACY  OF  THE  PRESSURE-ENERGY  ITERATION 

Consider  two  simultaneous  equations  in  two  unknowns 

y  *  /<*) 

x  «  g(y). 

A  procedure  is  to  be  devised  to  solve  the  equations  when  the  functional  form  of  /  or 
g  is  not  known  a  priori.  An  iteration  may  be  defined  by 

xt  =  g[/(x0)] 
xa  “  gt/(xx)] 
xi+l  *  gC/(x1)3 

where  *o  is  a  first  guess  and  is  the  1th  approximation.  The  iteration  converges  if  the 
derivative  of  the  combined  function  g/is  less  than  one  in  absolute  value. 


I*#*!  -  uri  <  1.0 


(See  Courant,  Differential  and  Integral  Calculus,  Vol.  I ,  p.  3S8.)  If  the  derivative  is  greater 
than  I,  an  iteration  using  the  inverse  of  the  two  functions  will  converge. 

If  is  the  error  of  the  i1**  approximation,  the  error  of  the  next  approximation 

is 

**fi  m  «'/' A* i- 

In  PISCES  IDL  and  2I)L  the  updated  pressure  and  internal  energy  density  are  related 
by  a  pair  of  simultaneous  equations. 


E"*1  -  /<„** 


*>  -  -  j(pml  +  pn)AVn^  +  En 


!»•*  -„n+i  .nd. 

P  *  g(V  ,  E  ) 

where  the  superscript  is  the  cycle  number  and  g  is  the  equation  of  state.  Following  the 
notation  above, 

g'  *  r/v 


Preceding  page  Milk 
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r  av 

2  V 


The  Gruneisen  ratio  is  usually  near  2  or  smaller.  Therefore,  the  error  after  one  iteration 
is  about  -AV/V  times  the  error  of  the  first  guess  and  the  error  after  two  iterations  is 
about  (AV/V)2  times  the  error  of  the  first  guess. 
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Notation: 


APPENDIX  E 

RELATION  BETWEEN  €r  <t  AND  COMPRESSION 
IN  PLANE.  CYLINDRICAL  AND  SPHERICAL  GEOMETRY 


R  =  undo  formed  coordinate  of  material  point 

r(R)  =  coordinate  of  material  point  in  strained  state 

u{R)  s  r(R)-R  *  displacement  of  material  point 

e  *  du/dr-  radial  strain 

r 

C  =  compression 


e  =  tangential  strain 


Plane  Symmetry 


unstrained 


r(R  ♦  JR) 


strain*  ti 


AV  »  [final  volume]  -  [inicial  volume] 

-  [c(R  +  dR>  -  r(R>]  -  [(R+dR)  -  R] 
«  r(R)  +  r'<R)dR  -  r(R)  -  dR 

-  [r,(R)  -  l]dR 
«  u/(R)dR 


AV  u'fRjr’’’. 

c  .  .  - - iUL,  -  -u  (R)  -  -  er 


Gj,  ■  0. 
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Spherical  Symmetry 


£V  -  [final  volume]  -  [initial  volume] 

*  [r3  (K  +  dRl  -  r3  (R)  ]  -  [(R+dR?  -  R3  ] 
r3  (R+  dR)  -  [r(R)  +r/(R)dR]3  -  r3 (R)  +  3ra(R)r  (R)dR 
(R+dR]P  -  R3  +  3Ra JR 

£V  «  [r3 (R)  +  3ra(R)r'(R)dR  -  r3(R)]  -  [R3  +  3RsdR  -  R3  ] 


-  3ra(R)r'(R)dR  -  3RadR 
»  3[u(R)  +  R]3[u'(R)+l]dR  -  3RadR 

*  [3uS(R)u'(R)  +6u(R)Ru'(R)  +  3Rau'(R)  +  3ua(R)  +  6u(R)K+  3K5  -  3Ra]dR 


fiV  [3uau'  +  6uRu '  +  3Rgu '  +  3ua  +  buK]dR 
v  3R*dR 


f  /  L  u®u '  2uu '  u£  2ul 

T  +ir  +  —  +  F+  Rj 

r  Sr+1  er+1l 

.^r+2u_+u.-_  J 


cr+1  u2  er+  1 


C  *  -  er  -  2et  . 


<3»> 

(3b) 

(3c) 
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APPENDIX  F 


FORMULAS  FOR  SOUND  SPEED  IN  ELASTIC-PLASTIC  MATERIAL 


‘v<l.i!l>ili 

v  =>  'ipccific  volume  =  l/p 
vQ  »  initial  specific  volume  *  l/po 
V  *  v/v0  ■•  relative  volume 
M  “  (vQ  -  v)/v 
C  -  (v0  -  v)/vc 
PQ  *  initial  density 
p  »  true  density 

The  definition  for  the  longitudinal  adiabatic  sound  speed  is 
/ 


>  \dp  / 


»  V 

P  ~  v„ 


jaz  az  ac  av  az  /  1  \  7  va  \  az  1  /  v  \2  l(2  i  az 

5-S535-5^|  (' v5, 

az  „  az  dy  Sv  „  az  /  M/  va  \  a  _i_  az 

op  “  du  dv  ap  *  au  \  v2/  \  p 0v0; *  p0  a^i 

c  .  (£\hmJ±M$  m  /J.  iS\* 

*  \ap/  v\p0  ac/  \p0  ap/ 
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